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ABSTRACT 



Context. Over the past decade a new generation of chemical models, in addition to the gas, have included the dust in 
the treatment of the interstellar medium (ISM) . This major accomplishment has been spurred by the growing amounts 
of data on the highly obscured high-z Universe and the intriguing local properties of the Solar Neighbourhood (SoNe) 
of the Milky Way (MW) Disk. 

Aims. We present here a new model able to simulate the formation and evolution of dust in the ISM of the MW. The 
model follows the evolution of 16 elemental species, with particular attention to those that are simultaneously present 
in form of gas and dust, e.g. C, N, O, Mg, Si, S, Ca and Fe. In this study we focus on the SoNe and the MW Disk as a 
whole which are considered as laboratories to test the physical ingredients governing the dust evolution. 
Methods. The MW is described as a set of concentric rings of which we follow the time evolution of gas and dust. 
Infall of primordial gas, birth and death of stars, radial flows of matter between contiguous shells, presence of a central 
bar, star-dust emission by SNae and AGB stars, dust destruction and accretion are taken into account. The model 
reproduces the local depletion of the elements in the gas, and simultaneously satisfies other constraints obtained from 
the observations. 

Results. The evolution of the element abundances in the gas and dust has been well reproduced for plausible choices of 
the parameters. The Mg/Si ratio, in particular, drives the formation of silicates. We show that for most of the evolution 
of the MW, the main process for dust enrichment is the accretion in the cold regions of the ISM. SNae dominate in 
the early phases of the evolution. We have also examined the main factors controlling the temporal window in which 
SNae govern the dust budget both in low and high star forming environments. The role played by AGB stars is also 
discussed. We find that IMFs with regular slope in the range of massive stars better reproduce the observed depletions. 
Conclusions. The classical chemical models nicely reproduce the abundances, depletion factors and dust properties of 
the SoNe and the main ingredients of the models are tested against observational data. The results obtained for the 
SoNe lead us to safely extend the model to the whole Galactic Disk or galaxies of different morphological types. 

Key words. Galaxies - Dust; Galaxies - Spirals; Galaxies - Milky Way 



1. Introduction 

In the fascinating subject of the origin and evolution of 
galaxies, the interstellar dust is acquiring a primary role 
because of its growing importance in the o bservations of 



the high-z Univers e (Omon t et al.l 12 001: Shaple v et al. 



200lUBertoldi et al.ir2003: Robson et al. 2004: Wang et al. 



2008albl : iGallerani et al.l 120101: iMichalowski et al.ll2010bllar ) 
and the theoretical spectro-photo metric, dynamical 
and chemical modelin g of galaxies ( Schurer et al. 200' 



Narayanan et al. 2010; Jonsson et al 



2010HPipino et al...2011:.PoDescu et al.l 



201fl[ 
20Tlfl . 



20101: Grassi et al 



strongly affects the production of molecular hydrogen and 
the local amount of UV radiation in galaxies thus playing 
a strong role in the star formation process via the cooling 
mechanisms (jYamasawa et al.ll2011|) . The inclusion of dust 
in the models leads to a growing complexity and typically 
to a much larger set of parameters influencing the results 
of the simulations to be then compared with the observa- 
tions. Indeed several question must be addressed, each one 
easily expanding the model: who are the main Stardust 
injectors in the interstel l ar medium (ISM) (jGail et al. 



Indeed, the evidence of highly obscured QSOs and 
galaxies already in place at high-z leads necessarily to a 
new generation of theoretical models where dust is a key 
ingredient that cannot be neglected, if we want to obtain 
precious clues on the fundamental question about when 
and how galaxies formed and evolved. First of all, dust 
absorbs the stellar radiation and re-emits it in the infrared 
deeply changing the shape of the observed spec tral energy 
distributions (SEDs) of obscured galaxies (jSilva et al.l 
Il998t iPiovan et all l2006t iPopescu et alll201ll) : second, it 



20091: IVahante et al.l 12001 iGaU et al1l2011at IPiovan et al. 



2011al )? Ho w much dust do they produce and on which 



timescales (jDraind [2009t iDwek et al.ll2009D ? What is the 
contribution and the role played by the molecular clouds 
(MCs)-gr own dust that form in the cold dense reg ions of 

1? How 



(IVHJs)-gr own dust that torm m the cold dense reg ie 
the ISM (iZhukovska et al.l 120081: IPipino et al.l 120111)? 



much dust is destroyed by SN ae shoc ks (jN ozawa et al."2006i 
200l iBianchi fc Schneided [20071: Ijones & Nuth 2011)? 
What is the typical minimal set of dust grains whose 
evolution should be followed and what could be a minimal 
set of dust grains to be used for a satisfactory description 
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of the chemical or spcctrophotonietric properties of the 
galaxy? 

To answer all these questions the theoretical models 
must include (i) a set of grains with suitable composition 
and properties and/or an ISM made of gas and dust in 
which the abundances of the elements are followed, (ii) 
a recipe for their formation/accretion and destruction 
in the ISM and (iii) a prescription for the yields of dust 
by the stellar sources. The duty cycle of the dust ca n 
be schematically summarized as follows (jJonesI l2004t) . 
Stars, mainly AGBs and SNae, inject material in the ISM, 
mainly in form of gas, but with a variable amount that 
condenses into the so-called star-dust. Once injected into 
the ISM, star-dust grains are subjected to destruction 
processes that restitute the material to the gaseous phase. 
However in cold and dense regions dust can accrete on the 
so-called seeds: the competition between the accretion and 
destruction processes, mainly via shocks, determines the 
total budget of dust in the ISM and the observed depletion 
of the elements that are involved in the formation of dust 
grains (jDwek 1998). Dust accretion mainly occurs in the 
very cold molecular clouds (MCs), where it induces strong 
cooling thus leading to the formation of new stars. The 
stellar winds from AGB stars SNa explosions more and 
more enrich the ISM with new metals and star-dust grains 
that are able to survive to the local shocks caused by SNa 
explosions. 

The Milky Way (MW) is th e ideal laboratory to ou r dis- 
posal to study the dust cycle (jZhukovska et al.l 120081 ) and 
its impact on the wider subject of galaxy formation. For 
obvious reasons, the MW provides plenty of observational 
data to which we can compare theoretical predictions, 
thus setting useful constraints on theoretical simulations 
and highlighting the role of the most important physical 
quantities leading the whole problem. Once this important 
step is accomplished, our modelling of the role played by 
dust can be extended to other galaxies such as local disk 
and spheroidal galaxies, high-z galaxies and QSOs. 

Starting from these considerations, in this paper we 
simulate the formation, evolution and composition of dust 
in the MW, both locally in the Solar Neighbourhood 
(SoNe) and radially along the Galactic Disk. We build up a 
detailed chemical model (the theoretical simulation is still 
the main tool to investigate the formation and evolution 
of dust in gala xies) starting from the pioneering study 
by iDwekl (1 19981) and takiii g into a c count the mor e recen t 
ones bv IZhukoyska et"ap (I2008D. ICalura et all (I2008D 
IValiantc et al.' ('2009'), 'Gall et al.' ('2011a'), 'Gall et alj 
(2011b), Ma t tsson (2011 ), Valiantc et al. ( 2011 
iKemper et all (|20nl ) and iDwek fc Cherchnefa (|20Tll 
The theoretical model w e are bu ilding up stems from the 
basic one with infall by IChiosil ([1980), however updated 
to the more recent version with radial flows of matter and 
presen ce of a central bar developed by iPortinari fc Chiosil 
(J2000t) . The stellar yields of che mical elements i n form 
of gas are those calculated by IPortinari et al.l ()1998D . 
The model follows the evolution of the abundances of a 
number of elements composing the ISM gas, includes the 
formation/destruction and evolution of dust and, finally, 
follows in detail also the abundances of those elements that 
are embedded in the dust grains. To this aims, the model 
makes use of the best prescriptions available in literature 
concerning dust accretion, destruction and condensation in 
the AGB winds and SNa explosions. These prescriptions 



have been already presented by I Piovan et al.l (|2011al ) to 
whom the reader should refer and will also be discussed in 
some detail here. 

The main test for any model of dust formation is 
given by the data on element depletion provided by SoNe 
of the MW to wh ich we will compar e our results. In a 
forthcoming paper (jPiovan et al.ll2011a ) we will investigate 
the radial dependence of chemical abundances and dust 
depletion across the Disk of the MW. 

The plan of the paper is as follows. In Sect. [5] we 
introduce the formalism and basic equations governing the 
temporal evolution of the gas, stars, and dust in an open 
model with radial flows of gas and dust for the galactic 
Disk of the MW. In Sect. |3] we summarize the current 
prescriptions for the star formation rate and initial mass 
function. In Sect. E] we introduce and describe in some 
detail the various processes responsible for the formation 
and growth of dust grains in the ISM, whereas in Sect. [5] we 
present the accretion rates into dust grains for a number of 
important elements we have considered. The yields of dust 
fro m AGB stars and SN a explosions are adopted according 
to iPiovan et al.l (J2011al) to whom the reader should refer 
for all the details. In Sect. JHIthe problem of dust destruc- 
tion by SNa shocks is faced. Then, in Sect. |7| we present 
the observational data for the elemental abundances in 
the Solar Neighborhood and define the reference set of 
abundances we have adopted. In Sect. |S] after summarizing 
the main ones between the many available parameters, we 
discuss and compare the effect of them on the formation 
and evolution of dust, at varying them between the many 
possible choices. In particular we examine the influence on 
dust of a different CO fraction in the ISM (Sect. JHH]), the 
effect of the IMF (Sect. HJ) and of the SF law (Sect. |531), 
the choice between different models for the accretion of 
dust in cold regions (Sect. 18. 4[) and, finally we discuss some 
interesting parameters however not discussed in our model. 
In Sect. [9] we present the results for our models of the Solar 
Vicinity in presence of dust and as a function of three 
important ingredients, namely, the initial mass function, 
the efficiency of star formation, the accretion time scale of 
primordial gas onto the system mimicking the evolution of 
the MW Disk. The effect of radial flows and central bar are 
always included according to the prescription developed 
in previous studies ( Portinari & Chiosi 2000) and also 
adopted in iPiovan et al.l (|2011cl) . The simulations are 
compared with the depletion of the elements in the SoNe 
under the constraints that we want that the local chemical 
properties are satisfied, such as the time evolution of the 
elemental abundances, the metallicity and iron enrichment. 
In Sect. (TUj we discuss the results we have obtained and 
draw some general conclusions. 

2. Chemical Evolution Model 

In classical models of chemical evolution, the Disk of the 
MW is subdivided in N concentric circular rings of a cer- 
tain thickness Ar, where r is the galacto-centric distance, 
in the case of plane geometry or N concentric cylindrical 
shells if the third dimension is considered. Each ring or 

shell is identified by the mid radius r^ with k = 1, ,N. 

In most cases, radial flows of interstellar gas and dust are 
neglected, so that each ring / shell evolves independently 
from the others. The physical quantity used to describe the 
Disk is the surface mass density as a function of the radial 
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coordinate r and time t: a{rk,t) is the mass surface den- 
sity at radius rfc and time t. Depending on the model, a can 
refer to the ISM (cr^), in turn split into dust or gas (cr^ 
or a*^ respectively), to the stars (a*) or to the total mass 
(simply a). At every radius r^, the surface mass density is 
supposed to slowly grow by infall of either primordial or 
already enriched gas and to fetch at the present age to the 
mass density profile across the Galactic Disk for which an 
exponential profile is best suited to represent the surface 
mass density distribution: a (rk,tG) ex exp {—rk/rd), where 
rd is the scale radius of the Galactic Disk, that is typically 
estimated of the order of 4 — 5 kpc. Since the final density 
profile is a priori known, one may normalize to it the cur- 
rent total surface mass density of the ISM " Ai" (sum of 
gas and dust), 

^,.,,0 = ^^. (1) 

CT(rfc,iG) 
Introducing the fractionary mass of the generic i-th element 
we have: 



M^{rk,t) = 



<^r{rk.t) 



= X^{rk,t)M{rk,t) (2) 



<7irk,tG) 

and therefore the fractional mass abundance Xii'''k,t) - 
Mi{rk,t)/Mirk,t), with Y.zXi{i"k,t) = 1. 

Similar expressions can be derived for the dust (indi- 
cated by " D) and the gas (indicated by " G" ) 

^^(rfe,t) 



Dirk,t) 
G{ru,t) 



(^{rk,tG) 
(T^{rk,t) 



(3) 

(4) 



crirk.tc) 

with (T^{rk,t) = a^ {rk,t) + a^ {rk,t) and a{rk,t) = 
0"^ {fkjt) + a* {rk,t), where a* {rk,t) is the surface mass 
density of stars. For single chemical elements we may write: 



D^{rk,t) 



X? {Trk,t)(TD{rk,t) 
^(rk,tG) 



G^{rk,t)^ 



Xi irk,t)aGirk,t) 



Xr{rk,t)D{rk,t) (5) 



^X?irk,t)G{rk,t) (6) 



(^{rk,tG) 

with J2i [xF {'''kit) + xf {^k^i)\ — 1, from which it follows 
that Y.^ xF irk,t) ^ 1 and E^ X? {rk,t) ^ 1. 

The fundamental equation describing the evolution of 
the ISM in abs ence of radial flows of ma tter between con- 
tiguous shells (jPortinari fc Ch iosi 2000) and processes of 
dust accretion/destruction (|Dwek.l998.) is: 



—Mi{rk,t) = 
at 



~Xi{rk,t)^p{rk,t) 

^J{rk,t-TM)R^{M)^{M)dM 



M, 



^M.(rK,i) 



(7) 



inf 



where </> (M) is the IMF and Mi and M^ are the lower 
and upper limits for the stellar masses, ip {rk,t — tm) is the 
star formation rate (SFR) at the radius r ^ and at the time 
t' = t-TM, Ri (M) = E.m/M (iPortinari et al. 1998) is the 
fraction of a star of initial mass M ejected back in form 
of the chemical species i-th. The three terms at the r.h.s. 



represent the depletion of the ISM due to star formation, 
its increase by stellar ejecta, and the increase by infall of 
external gas (either primordial or already enriched). 

Adding supernovae and radial flows. Type la supernovae 
originate in binary systems and have a fundamental role, 
in particular concerning the iron enrichment. The super- 
no vae rate and the ado pted formalism are the ones of 
Greggio fc Renzinil ()1983i ) and the formulation of equation 
des cribing chemical evolu t ion of the ISM is modified follow- 
ing [Matteuccr^^reggia (|i98^. The contribution of single 



stars, corresponding to a fraction (1 — A) of the total, is 
separated from the contribution of binary system, a frac- 
tion A of the total. Inserting the contribution of type la 
supernovae and integrating in time, instead that in mass, 
the equation for the evolution of the i-th component of the 
ISM is, 



Mi {rk,t) = -x»V'+ 



4> (M) R^ 



+ / ■ ^ 

/o 



+ {1-A) ^ 

Jt-TM„ , 



dM 

d-TM 



dt' + 



M(t) 



(|)R^ 






*«^|-0 



dM \ 

dTM J 

dt'+ 



dt'+ 



M(t) 



M(t) 



+ A 



*-'^Mi 



/(Afl) R^,l 



dMi 
' dr 



Ml 



dt'+ 



Afi(T) 



-I- RSNI ■ EsNI,i + 



+ 



jM.{r,,t) 



Mi{rk,t) 



inf 



rf 



jMdr.,t) 



(8) 



where V = '>Pirk,t'), (}> = 0(M), x = Xi(^fc,0) Ri = 
Ri (M), i?,,i = R, (Ml) and M {t - t') = M{t). The first 
term at the r.h.s. is as usual the one describing the de- 
pletion of interstellar material because of the process of 
star formation and it depends from the star formation rate 
and from the abundance of the i-th element considered. 
The next three terms represent the contribution of sin- 
gle stars to the enrichment of the i-th element. The fifth 
term is the contribution of the primary star in a binary 
system (assumed to be independent from the secondary 
star as far as it concerns the chemical yields). The sixth 
term is the contribution of type la supernovae. Finally, 
the last three terms are the infall rate of external gas, 
the outflow rate of matter due for example to the onset 
of galactic winds powered by supernovae explosions, and 
the radial flows of gas that determine the ISM exch ange 
between contiguous shells (jPortinari fc Chiosil l2000l ). re- 
spectively. Furthermore, /(Mi) is the distribution function 
of the mass of the primary star Mi in a binary system, 
between Mi^^in = Mb,i/2 and Mi^^ax = Mb,u, where 
Mb,i and Mb,u are the lower and upper limit of the bi- 
nary systems assumed respectively SM© and 12M0. Rsni 
is the rate of type la SNk and Esnia their ejecta of the 
i-th chemical species. M (r) = M {t — t') is the mass of 
a star of lifetime r, born at t' . It is worth noticing that 
various quantities depend on the metallicity Z (t) as well 
as on M: M (t) = M {t - t') = M {t - t' ,Z {t - t')) and 
Ri (M) = Ri (M, Z {t — t')) as stellar lifetimes and ejecta 



L. Piovan et al.: Formation and evolution of the dust in galaxies 



depend on metallicity. Ri{M) a re calculated on the ba se of 
the detailed stellar yields from iPortinari et all (|l998f ) and 



keep track of finite stellar lifetimes (no instantaneous recy- 
cling approximation). Eqns. ([5]) govern the evolution of the 
ISM. 

Separating gas from dust. For our purposes we need 
to formulate the equations governing the evolution of the 
dust in the ISM. Separating the ISM in gas and dust, the 
equations governing the evolution of the generic elemental 
species «-th in the dust are 



term describing the accretion of grain onto bigger particles 
in cold clouds; (4) the destruction term taking into account 
the effect of the shocks of SNae on grains, obviously giving 
a negative contribution. The infall term in the case of dust 
can be neglected because we can safely assume that the 
material entering the galaxy is made by gas only without a 
solid dust component mixed to it. 

Finally, from the equation for Di {rk,t) we can get the 
equation describing the evolution of the gaseous component 
Gt{rk,t), where Gi{rk,t) = Mi{rk,t) - D,{rk,t): 



^A(r„t) = 


-X 


°v+ 






I 


■ » / dM 


+(1-^) 


/ i' 


^5Z,R. 


J 


t-r, 


-fR.i 





dt'+ 

M(t) 

dM\ 

dTM J 



dt'+ 



M(t) 



+ (1-^) 



*-TA/p 



t-TM 



SNsc 
II 



+ A V 

■'i-TMi.min 

+ RsNlEsNI,iSc,i + 

' d 



^ 



Rr 



:,i-I^r ' 

dM 

dTM 



dM\ 

dTM J 



dt'+ 



M(t) 



dt'+ 



M(t) 



fiM,)6ZR.^,-(-^ 

dTMi 

dMi 



dt'+ 



M(t) 



/(Mi)5™,7?,,i 



drn 



dt' + 



M(t) 



dt 



Di{rk,t) 



+ 



D^{rk,t) 



+ 



rf 



+ 



^A(..,t) 



^A(r„.) 



(9) 



SN 



where xF — xF ii^k,t)- The first term at the r.h.s. of Eqn. 
([9]) is the depletion of dust because of star formation that 
consumes both gas and dust (uniformly mixed in the ISM) . 
The second term is the contribution by stellar winds from 
low mass stars to t he enr i chmen t of the i-th component of 
the dust. Following Id wekl ()1998f ). we introduce the so-called 
condensation coefficients Sf^ that determines the fraction 
of material in stellar winds that goes into dust with respect 
to that in gas (local condensation). The third term is the 
contribution by stars not belonging to binary systems and 
not going into type II SNas (the same coefficients S'^^ are 
used). The fourth term is the contribution by stars not 
belonging to binary systems, but going into type II SNae. 
For the condensation efficiency in the ejecta of type II SNas 
we introduce the coefficients S^\, the analog of 5™j. The 
p ossible choices for thes e coefficients are discussed in detail 
in iPiovan et al.l ([20113) ■ The fifth term is the contribution 
of massive stars going into type II SNae. The sixth and 
seventh term represent the contribution by the primary star 
of a binary system, distinguishing between those becoming 
type II SNae from those failing this stage and using in each 
situation the correct coefficients. The eighth term is the 
contribution of type la SNae, where again we introduced the 
condensation coefficients S^.^ to describe the mass fraction 
of the ejecta going into dust. The last four terms describe: 
(1) the outflow of dust due to galactic winds (in the case 
of disk galaxies this term can be set to zero); (2) the radial 
flows of matter between contiguous shells; (3) the accretion 






-XG,»V' + 



i~^Mq 



t-TMf, 



(1 - 5Z) R^ 



+ / V' 

Jo 





0(1- 


' sZ) R^ 


"SAfa, '- 




0(1-5 


^i) Ri 


/ dM 

V dTM 



dM 

dTM 

dM \ 

dTM J 



dt'+ 



M{t) 



dt' + 



M(t) 



dt'+ 



M(t) 



+ 






^ 



0(Af)(l-C.)^-(-^ 



dt'+ 



A 



A 



^ 



t-TAfs„^ 



/(Ml) (1 - 5,",) R,, 



^ 



f t — TMi ,Tni: 

+ RsNlEsNI,i ( 1 — S^i ) 4- 



/(Ml) (1-5™,)^,,! 

d 



dt 



G^{rk,t) 



M{t) 

dMi 

dTMi 

dMi 

dTMi 

+ 



dt'+ 

M{r) 

dt'+ 



M(t) 



inf 



^G.(..,t) 



D,{rk,t) 



+ 



+ 



^G. (..,.) 



+ 



rf 



D^(rk,t) 



(10) 



where again the outflow term — [Jt^i (''fci^)] '^^iH be 
fixed to zero because we do not have galactic wind for 
spirals with continuous star formation. Since the primor- 
dial material is likely dust-free we have [jjA^i (^7^)] ■ — 

It is worth noticing the following point: the stellar mod- 
els, upon which are based our yields, predict that stars 
with mass higher than QMq go into SNs, whereas those 
with mass lower than 6Mq first become AGB stars and 
later White Dwarfs. We must therefore split the third and 
the fifth member of Eqn. (|S]) in two parts, both in Eqns. 
(fTO)) and (IH), because the minimum mass dividing the in- 
tervals of AGB and/or SN^ belongs to the mass interval 
(3 — I2M0) describing binary systems going into type la 
SNaei. 



^ For example in Eqn. ([9|, using Scj. to indicate the generic 
condensation coefficient we have the following split: 



{I -A) 



is divided into: 



(1-A) 



t-TA/p 



t — T"AJq 



^ 



„ ,,,, / dM\ 



dt' 



M(r) 



^ 



t — TM 



B,l 



Sc,iRi 



A) I i, 



iRi 



dM \ 
dTM ) 

dM \ 

dTM J 



dt' 



M(t) 



dt' 



M(t) 
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To summarize. Indicating the contribution to the yields 
by stellar winds and type la and II SNae with the symbols 
Wi^Dirk,t), Wi^G{rk,t) and Wi^M {rk,t) (they can easily 
be reconstructed by comparison) and neglecting the outflow 
term, Eqns. dH]), ^ and dTU]) become: 



dt 



M^irk,t) = -Xi{rk,t)^p{rk,t) + +W,,Mirk,t) 



i^- 



{Tk.t) 



(11) 



rf 



3. The Star Formation Laws and Initial Mass 
Functions 

Star Formation. The law of Star Formation (SF) is a key 
ingredient of any model of galaxy formation and evolution. 
Unfortunately it is poorly known, so that many prescrip- 
tions for the SF rate can be found in literature. In this study 
we have consi dered several well known . SF laws adopted for 
the MW fsee lPortinari fc Chiosilll999l for details). 

A very popular prescription is the iSchmidtl (|1959() law. 
In our formalism it becomes: 



dt 



D^ irk,t) = -xf (rfe,t) V(rfe,t) + +W,,g {rk,t) 



+ 



J,D.drk,t) 

d „ , 
ja.irk,t) 



^A(r-fe,i) 



SN 



(12) 



rf 



dt 



Gi (rk ,t) = -XG.t (ffe ,t)'iP{rk,t) + W^.G {rk , t) 



^A(r„t) 



jG.irk.t) 



inf 



^^D.irk.t) 
^a.irk,t) 



SN 



(13) 



rf 



It is soon evident that the dust creation/destruction and 
the radial flows make the system of diffe rential equations 
more complicated than the original one bv lTalbot fc Arnefra 
(jl975|) for a one-zone closed-box model. As the ISM is given 
by the sum of gas and dust, only two of these equations 
are required, furthermore Eqn. (|TT|) can be used only if 
gas and dust flow with the same velocity. To proceed fur- 
ther we must now specify the law of star formation, the 
IMF, the stellar ejecta and the various rates describing gas 
infall, dust accretion/destruction, and radial flows/bar ef- 
fect. No details will be given about these ones. They are 
included into the model and they are mainly useful in or- 
der to reproduce the radial gra dients of abundance in the 
M W. The reader s hould refer to iPortinari fc Chiosil ()2000[ ) 
and lPiovanet al.l (|2011ct) . 



where Msnic is the separation mass that tell us if we must use 
the condensation coefficients of stellar winds (condensation of 
dust in the envelopes of AGB stars) in the mass interval between 
Mb,i and Msnx or the condensation coefficients of supernovae 
between Msnx e Mb,u- In the same way: 



A 



splits itself into: 



f{Ml)Sc,^R^,l 



dMl 
dTMi 



dt' 



Mi(t) 



A 



' — T"Mq 



f(M,)Si!M.,i{-'^^" 



+ A 



t-TAfo 



t-T"A/j, 



/(Ml) 4:>fi».i 



dMl 

dTM^ 



dt' 



M(t) 



dt'. 



M(t) 



*(rfe,t) 



dG{rk,t) 



dt 



(^{rkyto) 



a{rQ,tG) 



G^r^,t) (14) 



where the normalizat ion factor is a(rp,, tp) ^'^ ^' and v 
is in [t~^]. Following IPortinari fc Chiosil (|l999l) we adopt 
K= 1.5. 

This simple dependence of the SFR can be complicated 
by including other physical effects. For instance, the SF 
suited to spiral galaxies such as the MW, may include the 
effect of gas compression by d ensity waves (Roberts 1969; 
IShu et al.lll972l: IWvse fc Si lk 1989; Prantzos fc Silk 1998fl 
or gravitational instabilities (Wang fc Silk. 1994 1 . We have: 



^(rfc,*) = iy 



r 



a{rk,tG) 



(^{rG,tG) 



G^{rk,t) (15) 



where v is always in \t ^] and k — 1 (jKennicuttl 119981 : 
IPortinari fc Chiosilll999ll . 

Another possibility is to describe the SF as a balance 
between cooling and heating processes, that is the grav- 
itational settling of the gas ont o the Disk and the energy 
injection from massive stars (jTalbot fc Arnetnll975l : [DoDital 
119851 : JDopita fc Rvdeijll99^ . In our formalism, we have 



'i'{rk,t) = 



\rk,t)a"'-\rk,tG) 



c^(rQ,tG)"^ 



G'"(rfe,t) 



(16) 



where n = 1/3, m = 5/3 (jPortinari fc Chiosil |1999() and 
V i n \t~^]. This formulatio n is similar to the original one 



by Talbot fc Arnctt (1973) thus leading to similar results 
(IPortinari et al. 1998) . 

Initial Mass Function. The initial mass function (IMF) is 
perhaps the mos t imp ortant ingredient of chemical models 
of any kind (see iKrou pa 2002b, for a recent review of the 
subject). Brown dwarfs and very low mass stars whose life- 
times are longer than the age of the Universe, in practice 
lock up forever the chemical elements present in the ISM at 
the age of their birth, whereas intermediate and high mass 
stars of short life continuously enrich the environment with 
the products of thermonuclear reactions, thus driving the 
chemical evolution of the host system. They are also the 
factories of star-dust to be injected into the ISM by SNae 
explosions and strong stellar winds. The adoption of an IMF 
has two effects worth being mentioned here. First of all, a 
different slope of the IMF in the intermediate-high mass 
range would imply a different relative population of the 
stars contributing to the dust yields. Second, the net yield 
of metals and dust per stellar generation varies. According 
to its definition (see fo r instance iTinsley 1980 ; Pagcl 199J ; 
IPortinari et al.l l2004al) the net yield is the the amount of 
metals globally produced by a stellar generation over the 
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the fraction of mass locked up in living stars and remnants. 
Therefore, efficiency of metal and dust enrichment depends 
not just on the amount of metals produced per unit mass 
involved in star formation but on the ratio between this and 
the mass that remains locked in remnants or ever-lived low 
mass stars. The locked-up fraction, is therefore as crucial to 
the metal and dust enrichment as is the absolute number 
of the high-mass stars directly responsible for the produc- 
tion of dust itself. In a given model of fixed total mass, it 
is clear that an IMF bending down steeply at low masses 
will lead to a different locked up fraction with respect to 
a power-law, low-mass oriented IMF. The metal and dust 
production is accordingly affected. 

In this section we shortly present the IMFs we have in- 
cluded in our model of the Galactic Disk and SoNe. For the 
purposes of our study the IMF is assumed to be constant 
in time and space. All the IMFs are normalized assum- 
ing that the total mass encompassed by the IMF from the 
lower. Mi, to the upper, Mjj, mass limit of star s is equ al to 
1 M(7i. To this aim, following iTalbot fc Arnet"^ (|1975D and 
iBressan et al.l ()1994[ ). we define the parameter C, which de- 
scribes the fraction of total mass in form of stars stored in 
the IMF above a given mass M*. In other words, M* is the 
minimum mass contributing to the nucleo-synthetic enrich- 
ment of the ISM over a timescale of the order of the galaxy 
life 



nAf, 



c = 



j;;y{M)dM 



nAf, 



j::y{M)dM 



(17) 



This equation, at varying ^ and for fixed Mu and M*, 
can be reversed numerically to determine the lower limit of 
the distribution M^. The IMFs included in our model are: 

- The Sal peter IMF. Salpeter-like IMFs are very popular. 
These are an extensio n over the desi red mass range of the 
original Salpeter IMF (jSalpeteiJIlQSSD . This IMF is (M) = 
CsM^^-^^ with Cs de pending on the value of C- For a mass 
range [0.1 - 100] M© ((Portinari et al.ll2004al) we have Cs = 
0.1716 and a C = 0.3925. 

- The Kroupa IMF. In a series of papers Kroupa revised 
and updated the power -law IMF with a set of continuous 
multi- sl ope po wer-laws (jKroupa et al.l[l993t lKrouDall200ll 
I2002bl lal. l2007l e.g.). In the following we c onsider two cases. 
First, the IMF derived by Kroupa (1998) for field stars in 
the SoNe and used by [Portinari et al. (2004b). This IMF 
is typical of models of chemical evolution of disk ga.laxies 
(iBoissier fc Prantzoslll99l 120001: [Prantzos fc Boissieij|2000t 

iHou et al.ll2008h ^ . For Ml = O.lM p and Mu = l OOMr^, 
we obtain ( = 0.405. Second, the iKroupal (|2007D IMF, 
where taking Ml = O.OIMq and Mu = 100M(^, we get 
C = 0.38, slightly lower than in the above iKroupal (|1998[ ) 
because of the lower li mit extended t o brown dwarf regime. 

- The Larson IMF. iLarsonl (|1998D proposed an IMF in 



which the relative percentage of very low mass stars and 
sub-stellar objects is decreased due to the presence of an 
exponential cut-off. As a consequence of this there is a neg- 
ligible contribution to the locked-up mass, and in contrast 



^ Along this line it worth recalling that IMF with slope 
(1.6 ~ 1.7) in the high mass range , i.e. steeper than the 
Salpeter value, is from IScalol (|1986l ) and it is widely used 
in literature in chemical models of the MW even with dust 
llMatteucci fc Francqidll989l: iChiappini et al.|[l997l: lDweklll998l : 



a very high net yield per stellar g eneration, and a high pro- 
duction of metals and dust. The ILarsonl (|l986t ) IMF is 



,(A/) = CiM"i-35exp --^ 



Ml 
M 



(18) 



In practice this IMF recovers the Salpeter IMF at high 
masses, whereas at low masses the exponential cut-off de- 
termines the steep downfall after the peak mass Mp = 
Ml/1.35 = 0.25Mo. For a typical mass range [0.01 - 
IOOJMq we have C = 0.653 thus allowing for a high num- 
ber of intermediate-high mass stars. For the present aims, 
we will simply keep Ml constant with time or metallicity. 
We also consider the possibility for a modified Larson IMF 
adapted to the SoNe, in which the slope in the power-law 
factor for the h i gh m ass range has the value Af ^^•^, ac- 



cording to IScalol (jl986t ) and in agreement with recent IMFs 
proposed by Kroupa (see above). In this last case with the 
same mass range [0.01 — 100], we have C = 0.5, lower than 
the other case. 

T he Cha brie r IMF. Along the line of thought of lLarsonI 
, IChabried (|200lD proposes: 



(f){M) = CcM 



-2.3 



exp 



Mc 
' M 



1/4' 



(19) 



The two parameters in equation (|19l) are tuned on local field 
low mass stars and the functional form is proved to be o 
be valid down to the brown dwarfs regime (Chabrier 20021 ) • 
With a mass range [0.01 - 1OO]M0 and Mc = 716.4 we get 
Cc = 40.33 and ( = 0.545, not as high as in the Larson 
IMF, but still leading to high net yield and low locked up 
mass fraction. 

- The Kennicutt IMF. The iKennicuttI (|1983D IMF is used 
in literature to describe the global properties of spiral galax- 
ies and it is inspired by the observations of Ha lumin osities 
and equivalent width in external galaxies JK ennicutt et al.l 
11994 iPortinari et al.|[2004bHSommer-Larsedll996[) . With a 
mass range [0.1 — 100] M©, we get ( = 0.59. 

- The Arimoto IMF. This t op-hea vy IMF has been sug- 
gested by lArimoto &: Yoshiil (|1987f) to simulate elliptical 
galaxies and is introduced just for the sake of comparison as 
an extreme case: 4>{M) = CAriM~^'^. With a mass range 
[0.1 - 1OO]M0, we get C = 0.5. 



4. Growth dust grains in the ISM 

Dust grains form by a number of physical processes (see 
below) and once in the ISM they may grow in mass by 
accreting a number of atoms or molecules. In the following 
we examine in some detail the accretion mechanisms for 
which two different models are proposed. 



4.1. The Dwek (1998) model or case A 

iDwekl (|1998D view of grain growth can be summarized as 
follows. Let us consider an ISM whose dust component is 
formed by grains made of a single element z-th with mass 
TOgr,i • Let rrii be the mass of an atom (or molecule) of the el- 
ement z-th in the gaseous phase. Ni is the number of atoms 
(molecules) of the element z-th locked up into the mono- 
composition grains that we are considering and ngr^i the 



llMatteucci &: i'tanmiaiiaSiJI: lUniappmi et al.llliJbli'l: lUweidli^iJai : cumpoBuiun gianis unau we are cunsiueimg anu iigr^i me 
iRomano et al.ll200(]l : iFrancois et al.ll2004l : ICalura et al.ll2008l ') density of grains of z-th type in the ISM, considered in this 
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model as a single phase, without distinction between dif- 
fuse ISM and cold molecular clouds. Let now Ui be the 
sticking coefficient, telling us the probability that an atom 
(molecule) of the element i in the gaseous phase binds to 
the grains increasing the number of atoms/molecules on it. 
Finally ngas.i is the number density the gas made of atoms 

f SKBT \^ . ^, 

IS the mean 

V Trmi J 

thermal velocity of the particles of the element i in the 

gaseous phase with respect to the dust grains. We get: 

dNi 



(molecules) of type i-th and v 



dt 



(XilKX IT'gr.i^^gaSji^ • 



(20) 



Multiplying both members of the equation by the constant 
mass rrii, multiplying and dividing the second member by 
the mass of one grain irigr^i, introducing the mean thermal 
velocity v and using the following relation ^ (|Dweld[l998l ): 

mingas,i = 2mHnH2 f 1 j^ j (21) 

where af (r, t) and Ci (r, t) are functions of position and 
time, we obtain: 



da. 



D 



^D 



1 



.D 



tM 



(22) 



dt Tifl 

For the dimensional consistency of the equation, the multi- 
plying quantity has the dimension of the inverse of a time 
indicated by t^^q 



1 



(23) 



Ai \ Trmi / 

It may be worth of interest to give an estimate of the in- 
volved timescale: by me ans of typical values of the involved 
quantities, Dwelg (Il998[ ) obtains ~ 3 x 10'* yrs. It is possible 
now to define the accretion time scale of our i-th element 
onto the grains as: 



1 



^M 



(24) 



that is the inverse of To.i multiplied by the fraction of the 
i-th element in the gaseous phase. Again, we drop the de- 
pendence from r and t. Dwek (1998) adopting the ratio 
<jf /(J-^ ~ 0-''' gsts Ti^accr ~ 10^ yrs which turns out to be 
significantly shorter than the lifetime of molecular clouds 
where accretion takes place. This lifetime is estimated to 
be of the order of tnc ~ 2 - 3 x 10^ - 10^ yrs. The life- 
time of a molecular cloud is comparable to the lifetime of 
the most massive stars born in it. These stars indeed in- 
jecting great amounts of energy by stellar winds and SNai 
explosions eventually disrupt the cloud. Similar estimates 
have been made in stud ies on the temporal evolution of 
dust in the ISM by Ca lura et al.l (pOOS). ~ 5 ■ 10^ yrs, 
and lZhukovska et al.1 (|2008l) . ~ 10^ yrs, and are consistent 
with observational a nd theoretical estimates of the lifetimes 
of molecular clo uds (|Matzneill20"ol iKrumholz et al.ll2006t 
iBlitz et al1l2007l see e.g.). One may argue that we are com- 
paring two different timescales, i.e. with ti^accr ^ tMC- 



The above relation says that the mass of the element i-th 
in the gaseous phase can be expressed by means of the density 
of the molecular clouds, i.e. the site in which the grains grow. 

2mHnH2 is the mass of the molecular cloud, ( 1 — (7^ /(Ji ) is 
the fraction of the element i-th in the gaseous phase as a function 
of the abundance of the same element in the dust and ISM. 



However, inside molecular clouds there are many physical 
processes continuously stripping atoms and and molecules 
from the grains, like UV radiation and cosmic rays. The net 
effect of it c ould be that Ti,accr is considerably lengthened 
(|Dweldll99a l. likely up to 6 • 10^ yrs, so that ti^accr — tmc- 
Furthermore, all the dust grains are in cold molecular 
clouds where the growing occurs. The simplest way to take 
this into account is to divide the above estimate for U^accr 
by the fraction of dust in molecular clouds. The resulting 
accretion timescale would be ti^accr ~ 1 — 2 x 10^ yrs. 

There is another important point to consider on which 
indeed several studies of the ISM evolution are based: both 
the destruction timescale (driven by SNs explosions) and 



the accretion timescale i,- 



have the same dependence 



on the star formation rate ip ^nd surface mass density ai. 
SNffi rate can be expressed as Ti^^nr oc a/ip. The lower the 
star formation rate, the lower is the number of supernovae 
explosions and finally the longer the destruction time of 
grains due to shocks. We clearly expec t Tj.gccr to v ary over 
the evolutionary history of the Galaxy: iDwe3 (|1998[ ) divides 
the accretion timescale by a^^'~^ /a and reasonably assumes 
that ip ex cr^^'^. In this way, Ti,accr oc atp. Since both Ti^accr 
and Ti_Knr have the same dependence from a /ip iDweld 



suggests that the ratio between the two timescales 
is constant: Ti^accr / Ti^snr = 1/2. If now we consider the ac- 
cretion/destruction terms in eqn. (J12p we have: 



dt 



(daf 



V dt 



.D 



^D 



n 



(25) 



The contribution of the creation/destruction terms is 
greatly simplified and the numerical solution of the dif- 
ferential equations describing the dust evolution becomes 
very simple. More refined is the solution adopted by 
ICalura et al.l (|2008[ ) . In brief, (i) only tq,; is fixed and n^accr 
remains an explicit function of the ratio a^ /ai\ (ii) no a pri- 
ori correlation between Ti^accr and Ti^snr is supposed to ex- 
ist, the two timescales are defined independently, and both 
plays a role in the af evolution. In our model we adopt the 
same strategy. 

4.2. The Zhukovska et al. (2008) model or case B 

A step forward in mod elling grain acc retion has been 
made by I Zhukovska et al.l (2008) . In brief. IZhukovska et al.l 
(|2008D 's picture can be summarized as follows: (i) grain ac- 
cretion is tightly related to molecular clouds and a multi- 
phase description of the ISM would be required; (ii) the 
evolution of chemical elements is calculated in presence of 
some typical dust compounds that are representative of the 
grain families growing inside the cold regions of the ISM or 
ejected by the stars, namely silicates, carbonaceous grains, 
iron grains, and silicon carbide. This model has two impor- 
tant advantages. First, it allows for different physical situa- 
tions each of which implying different timescales for each el- 
ements and different relationships between Ti^accr and Ti^snr 
which are independent; second, accretion in cold regions is 
described in a realistic way; third, the model can be eas- 
ily incorporated in a multi-phase description of the ISM to 
estimate the amounts of gas in the warm and cold phases. 



L. Piovan et al.: Formation and evolution of the dust in galaxies 



Let us now examine some of the key points and basic 
equations of the model and how they are adapted to our de- 
script ion, where contrary to what made in Zhukovsk a et alj 
(|2008l) we track the evolution abundances of single elements 
in dust instead of ad-hoc families of dust grains. 

The growth of the generic type of grains j-th is driven 
by the least abundant element indicated as element i-th 
(otherwise called the key element) among those forming 
the grain. All the other elements concurring to form the 
j-th compound adapt their abundances to that of the key 
element. For a mixture of accreting grains made of silicates 
(pyroxenes and olivines, i.e. magnesium- iron sulfates), car- 
bonaceous and iron grains, the key elements are Mg, Si or 
Fe for silicates, depending on which has the lowest abun- 
dance, C for carbon grains, and Fe for iron grains. Most 
of these elements accrete some specific atomic or molecular 
species called growth species. The growth in mass of a single 
j'-th type grain is: 

J o gw aD i,j gw /nc\ 

-^ = Sjay, A^mH^n] (26) 

where ruj is the mass of one grain of the j-th type, Sj is 
the surface area, aj is the sticking coefficient, v^^ is the 
relative velocity of the accreting species with respect to the 
grain, A^m-a is the mass of one atom (or molecule) that 
is accreting on the grain of type of type j-th, n^™ is the 
number density of the growth species j, and finally f^f T/^i^ 
is the ratio between the number of atoms of the key element 
i-th in the growth species for the dust grain of type j-th 
and the number of atoms of i-th element we need to build 
one grain of the j-th type. 

After some substitutions a nd simplifications for wh ich 
all the details can be found in IZhukovska et aD (|2008[ ). it 
comes out: 



dk 



dt 



Ji,j (.t Ji,j) 



where fij a p^ is the fraction of the key element i-th 
that is already condensed onto the dust type j-th, called 
degree of condensation. 1 — fi,j is therefore the fraction of 
the key element i-th still in the gaseous phase and available 
to form growth species for the j-th dust component. An 
explicit for mula for t^^ as a function of all the relevant 
quantities ([Zhukovska et al.ll2008l see Eqn. 31) is: 



ISM. If /o,i J is the initial degree of condensation of the key 
element i-th of the grains j-th and the molecular cloud will 
be destroyed after a lifetime t we have that the effective 
amount of material that is given back to the ISM after the 
cloud dissolution is: 



M. 



MC 



{t)^{f^At)~h^.J)xi 



M. 



MC 



(30) 



where /i.j (i) — /o.i,j is the fraction of key element condensed 
with respect to the initial condensation fraction, xfmax i^ 
the maximum possible fraction of dust of the kind j-th 
that can be formed in the molecular cloud and that would 
exist if all the material is condensed in dust of type j-th. 
Multiplying by the mass of the cloud we get the mass of 
newly formed material M^'^^ returned to the ISM. The 

fraction xf.max i^ given by 






A?e 



3 ""^.J 



(31) 



as shown in the footnote below ^ . 

Not all the MCs have the same lifetime: therefore there 
is a certain probability P (t) that a molecular cloud is 
destroyed in the time interval t - t + dt. According to 
IZhukovska et al.l (|2008[l the probability is well represented 
by an exponential law with time scale tmc- The instanta- 
neous condensation fraction at the time t of the key element 
i-th relative to the dust compound j-th is derived from inte- 
grating over time and considering for each cloud of lifetime 
X, born at i — a; the condensation fraction holding at the 
time of birth fo^ij (t — x). Since the lifetime distribution 
function quickly decays over time, only the MCs born some 
Tmc before t will contribute in practice. The instantaneous 
condensation fraction is 



(27) hAt) 



1 



fo,^.J{t~x)e^/^^e--/-"o 



TMC Jo l-/o,,j(i-x)(l-e^/^r) 



—dx. (32) 



This procedure is however time consuming when applied 
in practice. A much simpler approach is provided by using 

the average restitution mass i? and the average con- 
densation fraction /, ■ assuming fo.ij as a constant. Using 
Eqn. (1301) and integrating between and t — oo we obtain 
the mean value of fi,j: 



3 



:46- ^ '^ 



Af 






3.5 ■ 10" 



Myr . (28) J 



Integrating eqn. (1271) with the initial conditions /jj- — 
fo,i,j for t = we get the final equation describing the 
evolution of the condensed fraction: 

Jo,i,je ' 1 



k3 (t) = 



(29) 



1 - hr.3 + h.,3^'^^^ ' 

As in real conditions the grain accretion takes place in cold 
MCs with a finite lifetime tmc, the above accretion time 
scale must be compared with tmc- As already discussed for 
model A of accretion, the molecular lifetime constrains the 
real time interval during which the exchange of matter be- 
tween cold clouds and surrounding medium can happen. It 
is the interplay between the accretion/disruption timescales 
to determine if and how much material can be condensed 
inside the cloud before it dissolves into the surrounding 



1 



''' TMC Jo 1 - /o,.,, + /o,»j-e*/"^ 



lM22l^^^_e-'/^"''dt. (33) 



^ The number density nn of hydrogen atoms is given by 
nu ~ pMC I [(1 -I- 4eHe) ?tih] which in turn follows from the series 
of equalities: pmc ~ ninnu + Amuuuc = "Ih"-h + 47nH?iHeHo ~ 
mnnuil + 4:eHe), where ene = nuc/nn. The fraction X3,d,max 

_ P],,nax _ Aj mn{ti,jnYl) _ 



is given by Xj 



pMC 



■ mHnH(l + 4eHe 



— — — jy. The explanation is quite simple: dividing the 

(1 + 4,€He)Vr,j 

number of atoms eijriH of the key element available in the molec- 
ular cloud by the number of atoms v^j of the same key element 
that are tied up when we form one dust grain, we get the max- 
imum number of dust grains that we can form per unit volume. 
Knowing the mass of a single dust grain, Afnin, the the mass 
density pf^^nax immediately follows. 
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Inserting Eqn. p3p into Eqn. (15(11) we get AIj^mc, that is 
the time averaged mass of dust of type j-th returned to 
the ISM. However, this is only part of the story, because 
the molecular clouds have also different mass, so we need 
to evaluate their mean mass M mc ^^nd use this instead of 

Mmc- 

Fractionary mass of dust grains of j-type. With aid of 
Eqn. (j33p we can derive the fractionary mass of dust of type 
j-th produced by the accretion of grains in cold clouds per 
unit area and unit time: 



dt 



1 



tmc 



{.fij - fo,i,j) ■ Xj^max^MC ■ (34) 



Looking at Eqn. p4l) . it requires to know the amount 
of molecular gas of the galaxy we are going to model. 
To do this, a multi-phase description of the ISM taking 
into account the exchange of matter between cold and 
warm phases would be required. This is not possible with 
the present model because Eqns. ([5]), ©, and ([TU]) have 
been formulated for the single phase description and with- 
out a description of the gas exchange between cold and 
warm phases. To overcome this limitation some adjustment 
of the model are needed. Introducing xmc = o^^^ ja^ 
as the fraction of mol ecular gas and with some passages 
(jZhukovska et al.ll2008[ ) the final equation for grain growth 
in MCs in the frame of a single-phase ISM. 




Fig. 1. Absolute value of the differences between the /^ 
calculated integrating and averaging on the MCs lifetimes 
as in Eqn. (I33p and the same quantity der ived from the an- 
alytic al formula of Eqn. ((36|) suggested bv lZhukovska et all 
()2008D . A grid of 200x200 suitably discrete values of /o,»,j 



and tY bas been considered. 



(da, 



D 



Xmc 



(7 



D 
i,j j,max 



of) . (35) 



V ^* ) accT {^-XMc)tMC 

We need however to fix xmc '■ in IZhukovska et al.l ()2008[ ) is 
fixed to 0.2 in agreement to the observations for the SoNe 
at the current time. However this approximation could not 
probably hold for other phases of the evolution of the S oNe 
or other regions of the MW. In iPiovan et al.l ()2011ct ) , by 
means of the data available for the MW disk about molec- 
ular hydrogen, neutral hydrogen, total amount of gas and 
star formation we suitably extended xmc by means of 
Artificial Neural Networks (ANNs) in such a way to have 
the amount of MCs as a function of the current local vari- 
ables describing the system. T his recipe will be adop ted in 
this work. For more details see lPiovan et al.l (J201ld ). 

Degree of condensation. Another critical parameter^ to 
be examined with accuracy is the condensation degree /, .- . 
It depends on the comparison between tmC: the average 
lifetime of the MC, and the typical accretion time of grains 
Tj^ of Eqn. (pS)) . If Tj^ » tmc then we expect that only 
small amounts of dust are produced because the cold clouds 
will be destroyed before that dust has enough time to grow. 
On the contrary, if t^^ <^ tmc j we expect the condensation 
of dust grains to occur almost corn p letely and consequently 
dust to increase. IZhukovska et al.l ()2008 l) propose an ana- 
lytical expression for the condensation degree as a function 
of the ratio t^ /tmc (assuming a constant fo^i.j)'. 



■>i-J 



1 



/■2 

■'0,1,3 



i + ^W 



(36) 



According to IZhukovska et al.l (J2008|) , this expression is no 
longer valid only for very small values of /o,i j- • The situation 
for our models is more complicated because we are dealing 



with large variations of fo.i,j and t^^ . To clarify the issue, 
we checked in advance whether or nor the above relation can 
be still used. To this aim we have calculated the differences 
between f^j from t Eqn. ((36|) and the relation of Eqn. ((33|) . 
This comparison is made for t|"" < IOOOtmc and rj"^ > 
tmc/50. If Tj'' > IOOOtmc o r rf < tmc/50 we us e the 
Taylor expansion suggested bv lZhukovska et aTl (I2008D deal 
with the extreme cases of very slow or very fast accretion 
time. 

The results are presented in Fig. 14.21 for a 200 x 200 grid 
of models with uniform spacing of /o,ij and t^ growing 
linearly from low to high values. The error can be signifi- 
cant, in particular for very low /o,i.j and very short t^ or 
high t!-^ and average fo,ij- 

As we discussed in Sect. 14. 1[ t he sa me accretion term in 
iDwekl (|1998D and lCalura erall (|2008h is 



dai 



dt 



'^0,1 



Z^-^--k ■ (37) 



The basic differen ce b etween the two app roaches is that 
while iDwekl (|1998D and lCalura et al.l (|2008D with describes 
the e volution of the abundan ce of the element z-th in the 
dust, IZhukovska et al.l ()2008D follow the evolution of the 
abundance of the j'-th type dust in a given sample molecules 
considered to represent the situation of the ISM. To adapt 
the equations to our formulation for single elements (see 
Sect.H]), we start from: 

(38) 



n 



^M 



which summing up all the types of dust j-th in which a 
given element is locked up the element j-th is locked (in- 
dependently on whether or not it is a key element for that 
kind of molecule) becomes: 
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with 



M 



^D \~^ ,,D ^i.j „ 



D 



(39) 

(40) 



where as usual af is the surface mass density of the i-th 
elenient in dust, while a^ is the surface mass density of the 
j'-th type dust. The above relations says that the amount 
of element i-th contained in the dust is given by the sum 
of the contributions by all the j-th kinds of dust in which 
i is involved. For each one of the j'-th types, we divide af 

by the unit mass A^mn of one species j'-th and find the 
number density of dust of type j-th. Then, for each one 
of these molecules 7-th we have vP, atoms of the element 
i-th. Finally, we multiply by the atomic mass A^jmu of 
the element «-th to get the mass of the «-th element coming 
from the j'-th type dust. If we derive respect to time we 
obtain 



dt 



TOT 



E 



D ^i,j 

-■ Af 



'dal 
dt 



(41) 



TOT 



where the total evolution of dust is given by the sum of the 
variations due to the star formation, the Stardust injected 
in the ISM, the destruction of the grains by SNa shocks, 
and the accretion process in the cold regions. Taking into 
account only this last one: 



da, 



D 



dt 



ED ^i,j 
'■J A^ 



dal 
dt 



(42) 



The amount of element «-th contained in dust is given by 
the sum of the contributions by all the jth kinds of dust in 
which I is involved. For each j'-th type, we divide a^ by the 

unit mass Aj'rriH of one species j-th of dust and find the 
number density of the j-th type dust. For each j-th dust 
we have i^/' atoms of the element i-th. Finally, we multiply 
by the atomic mass Ai^mn of the element i-th to get the 
mass of this coming from the j-th type dust. Inserting the 
expression for (dcr^/di) we have: 



f da, 



D 



V dt 



u^ A- 



3 ^ 

Xmc 



dal 
dt 



aa 



3 



\D 



(1 -Xmc)tmc 



{Ik,]^j,max -<^j) ■ (43) 



Let's now introduce aj^^^^^, which is related to the amount 
of key-element fc-th in the dust of type j-th present in the 
ISM, i.e. a-^. We can express aj' by means of the key- 
element fc-th of the j-th dust that is locked up into the 
j-th component itself. This quantity is not known and not 
tracked by our model. All we know is how much of a given 
element is locked up into dust, but we do not know how 
much of that element is locked in every kind of dust. For 
example we are able to follow how much silicon is locked 
up into dust, but not how much silicon is stored in pyroxene, 
olivine and quartz. Doing the substitution we have: 



dot 
dt 



<3^^^^ 



Xmc 



Y ^f ^Mc (1 



Xmc) 



fk 



A^ 



\D 



k.j A D "k,j 



A, 



k,3'^k,j 



n 



(44) 



da. 



D 



dt 



Xmc 



X^ 



^1,3 



,D 



■cr (1 ~ Xmc) tmc ^ ^k.^v^ ^ 
• [Jk.j'^kj - <^k,j\ 
Posing K = Xmc/ {tmc (1 - Xmc)), we have: 



(45) 



daf 
dt 



K 



E(7 
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1,3 i, 3 



^,3 
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E 



A 



A 



kj'^kj 



{fk 






'Kj, 



(46) 



where k indicates the key element of the dust of type j-th, 
because it may occur that a dust molecule contains the el- 
ement i-th of which we tracking the evolution, but not the 
key element that could be of another type. For example, 
we could be following the evolution of carbon and consid- 
ering the silicon carbide that contains carbon, but not as 
key-element. As a consequence of this, the above has been 
split in two parts: in the first one we include the part of the 
j-th type dust grain in which the element i-th is locked as 
key-element] in the second part we have the dust molecules 
containing the i-th element but not as key element^ which 
instead corresponds to the element fc-th. Therefore the sum- 
mation goes for all the k ^ i. The notation a^ , is meant to 
indicate the surface mass density of the fc-th element locked 
up into the dust compounds of type j-th. 



5. Dust accretion rates for a few important 
elements 

At this stage we need to specify the dust accretion rate 
given by Eqn. (|46p for some specific elements whose evo- 
lution we intend to study. The model follows 16 elements 
(|Portinariet all [1998): H, ^He, ^^c, "C, "N, ^^N, i^q, 

"O, 180, 22Ne, 20Ne, 24Mg, ^SSi, 323^ 40^^^^ 56pg_ 

Four families of dust can be identified and tracked by 
the model: they are the silicates, carbonaceous and iron 
grains, and silicon carbide (it is worth clarifying that SiC 
does not form by accretion in MCs). We must distribute 
the 16 elements followed by the model in the four families 
of dust. First of all, H, "^He, ^''Ne, ^^Ne are noble gases 
not involved into the dust formation process by accretion. 
Second, since the hydrogen is by far the most abundant el- 
ement in the Universe, even if it may be contained in many 
dust molecules (such as for instance the PAHs), its abun- 
dance in the ISM gas will not be significantly affected by 
dust. For what concern ^ ^p^ i7q^ isq^ 2Ay^^^ 28 g^^ 56pg^ 

according to the simplified I Zhukovska et al.l ()2008[ ) scheme, 
they are involved in silicates (all of them) and iron grains 
(just some of them). For ^''N, ^^N, ^28 and '"'Ca there is no 
specific dust family to associate them. However, the obser- 
vations tell us that they are or could be depleted -see fo r 
instance the case of Calcium (|Whittetll200l lTielensll2005h - 
, so we must set a plausible yet simple scheme to describe 
them. Finally, ^^^C and ^'^C are both involved into the for- 
mation of carbonaceous grains and silicon carbide. In the 
following we will examine the case of ^^C, I'^C, "'^^O, ^^O, 
1^0, 24Mg, 28si and ^^Ye in some detail. 
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5.1. Carbon 

Carbon is present in dust as key constituent of many 
molecules containing carbon (C-molecules), otherwise 
named carbonaceous grains. It is also part of important 
dust molecules like silicon carbide (SiC) where, however, 
the key element ruling the process is silicon. SiC grains are 
mainly related to AGB carbon stars and for them there is 
no process of accretion in cold clouds of the ISM. Finally, 
the problem of the formation of SiC in the cjccta of super- 
no vae is still deba t ed. 

iNozawa et al.l (|2003[ ). studying dust formation in pop- 
ulation III supernova argues that they cannot form SiC 
grains with radii comparable to the pre-solar SiC grains 
identified as supernova condensates, because silicon and 
carbon are firstly locked into other carbonaceous grains or 
silicates before the formation of SiC grains can take place. 
However, we observe a small but significant percentage of 
X-type SiC grains (one of the minor types A, B, X, Y and 
Z, that is believed to be formed in the ejecta of type II 
SNae because of the isotopic signature) with respect to the 
mainstream SiC (about 90% of the total and whose ori- 
gin is attributed to AGB C-stars, because of their similar 
12(-i^13q ratio, signatures of s-processes and the 11.3 SiC 
/zm feature in C-stars) . The ratio between X-ty pe and main- 
strea m type can be quantified in a bout 0.01 (iHoDpe et al.l 



l2000tl . Therefore, as proposed by iDenea ult et al.l (|200 



special conditions and/or chemical pathways might have 
to be considered in order to realize the formation in type 
II SNae of SiC grains. To reproduce th e obse rved X-type 
over mainstream ratio, IZhukovska et al.l ()2008[ ) considers a 
small formation of SiC in Type II SN«, however this is just 
an ad-hoc recipe based on meteorites and does not throw 
light over the uncertain amount of dust, and X-type SiC, 
produced by SNae. 

The equation describing the accretion of carbonaceous 
grains of the ISM is: 



dal 
dt 



-^{/c,c(l-eco)<c-<c} (47) 



where ctq^q is the amount of the available carbon, ^co a 
factor accounting for the carbon embedded in molecules 
like CO that do not take part to the accretion process, 
(Jqq is the carbon already condensed into dust, and f q q 
is the average condensation factor (see Sect. U model B). 
The accretion of C-molecules in molecular clouds happens 
according to a time scale defined by Eqn. [55] which for the 
specific case of carbon is: 



46 



{Ain"^vS,c 






pE 



10-^ 



3.5-10" 



Myr. (48) 



This relation deserves some comments. First of all 1/^^ = 1 
because carbon is present in the cold regions of molecular 
clouds as single atoms or molecules with t ypically one 
atom of carbon. Second, p^ = 2.26 g cm~'^ ( Li &: Draind 
l200lh . To determine A^ we should average the masses 
of the various type of molecules present in the molecular 
clouds. Carbon is embedded in CO (with a percentage eas- 
ily up to 20 or 40%). There are also lots of C-molecules like: 
CO2, H2CO, HON, CS, CH4, C2H2, C2H6, HCOOH, OCS, .... 
We consider Aj3 = 28 corresponding to CO as a sort of 
mean value fairly representing the above group. A^ ~ 12 



assuming single atoms as the typical accreting species. 
Finally, we need to specify ec,c the abundance by number 
with respect to H of the key-species C driving the process 
of accretion of dusty C-molecules. We have: 



£c,c — (1 — Ceo) ■ 



^M 



.D 



12< 



(49) 



In this relation, the factor (1 — £,co) takes into account that 
fact that the carbon already in CO is not available to the 
accretion process. This correction is not negligible. Second, 
we subtract from the carbon atoms available for accretion, 
those ones already embedded in dust or belonging to un- 
reactive species like SiC. Finally, we obtain for the typical 
accretion time: 

^ Afyr. (50) 






1.80 



(1 -Cco)nH \a^ -a'(j 



'H 



A similar relation is assumed for ^'^C, in which the only dif- 
ference is the 2.03 instead of 1.80. This assumption stands 
on the notion that ^^C, much less abundant than ^^C, un- 
dergoes the same kind of accretion processes and forms 
similar C-molecules. Clearly this is a simplified view of the 
subject, however sustained by the fact that, on average, the 
mean N(-'^^CO)/N(-'^^CO) ratio in diffuse clouds is close or 
compatible with the local interstellar isotope ratio ^^C/^'^C, 
even if the N(^^CO)/N(^'^CO) can be easily enh anced or re- 
duced by a factor of two with a large spread ( Lisztl l2007t 
IVisser et al.l 120091 ). This allows to think that the relative 
percentage of carbon isotopes in molecules and dust should 
not be on average too different. 

To practically use Eqn. (|47|) . we need fo.c,c to get the 
average condensation factor fee- fo,c,c is the initial con- 
densation of the key-element C of the carbonaceous dust 
compound. If xf is the mass fraction of dust j-th contained 
in the part of the ISM that is not stored in cold clouds, we 
have: 



/o. 



c,c 



Xc 



Xc 



^C.max 



^C.max 



n 



(51) 



where to derive it, we simply replaced the mass fraction of 
dust in the ISM not belonging to cold regions with the mass 
fraction embedded in the cold ones: Xc ~ Xc- This simply 
means that after every cycle of MCs there is complete resti- 
tution and mixing with the ISM, so that on the average the 
ISM dust content always mirrors the initial mass-fraction 
already condensed at the formation of the MC. 



5.2. Silicon 

Silicon is present in many types of grains generically re- 
ferred to as silicates, one of the most complicated fami- 
lies of minerals. We consider only some kinds of silicates, 
commonly used in models of dust formation. First of all, 
we list the pyroxenes that are inosilicates (i.e. a family 
containing single chain and double chain silicates) in gen- 
eral indicated by XY(Si, Al)206, where X represents ions 
such as calcium, sodium, iron^+ and magnesium, whereas 
Y represents ions of smaller size such as chromium, alu- 
minum, iron^^, magnesium. Aluminum, while commonly 
replacing in other silicates, does not often do it in pyrox- 
enes. Typical pyroxenes used in dust models are the end 
members of the enstatite-hypersthene-ferrosilite series de- 
scribed by Mg^Fei-xSiOa with < x < 1 determining 
the partition between the two possibilities. Enstatite is the 
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magnesium end-member MgSi03 of the series, while fer- pyroxenes and ohvines in which silicates are subdivided, we 



rosilite is the iron end- member FeSiOa: both are found in 
iron and stony meteorites. Second, we have the olivines. 
They are a series of minerals falling in between two end- 
members, fayalite and forsterite. They can be described by 
[MgxFei_x]2Si04 with < x < 1 determining the partition. 
Fayalite is the iron-rich member with Fe2Si04. Forsterite is 
the magnesium-rich member with Mg2Si04. The two min- 
erals form a series where the iron and magnesium can be 
interchanged without much effect, also difficult to detect, 
on the crystal structure. Olivine can be found for example 
in many meteorites like the iron-nickel ones, in the tails of 
comets, and in the disks of dust around young stars. Both 
pyroxenes and olivines are expected to exist in the interstel- 
lar dust because they form from magnesium and iron that 
are abundant. In addition to this and limited to the calcu- 
lation of stellar dust yields, we included quartz, the most 
common mineral on the surface of the Earth, composed by 
Si02 and silicon carbide, SiC. Silicon carbide in expected 
to form in the ejecta of carbon-rich stars and type II SNae, 
but not by accretion in cold clouds. 

Two parameters control the pyroxene-olivine mixture: 
the Mg fraction x, assumed to be the same for pyroxenes 
an olivines, and the abundance ratio between the amount 
of magnesium and silicon bound in dust Af.j /Ag^. This 
latter determines the fraction of silicates with olivine sto- 



ichiometry Fqi — ^A/„/(^si ' 2;) — 1, while for pyroxenes 
we can define Fpyr = 1 — Fqi- The two parameters are 
fixed and chosen according to present day obs ervations of 
IR emission and abundances in the diffuse ISM (iDwek et al.l 
ll997HWhittetll20n3l:lD^^l2005l:IZhukovska et al.ll2008h . i.e. 
X = 0.8 and Af^g/Ag^ = 1.06 — 1.07. Varying these parame- 
ters does not significantly affect t he total e f ficien cy of dust 
production by the MCs (Zhukovska et al.l l2008h . For the 
purposes of this study, they can be kept fixed. The equa- 
tion describing the evolution of silicon is nearly the same 
no matter whether or not silicon is the key- element for py- 
roxenes and olivines. 

The Mg/Fe case. Let us first suppose that the key- 
element is magnesium/iron and not silicon. Using Eqn. (j46p 
we get: 



dig 
dt 
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(52) 



where fi^pyr and f^oi ^^^ the average condensation frac- 
tions of pyroxenes and olivines, obtained according to the 
procedure described in Sect. HI Furthermore, i^g p is the 
number of silicon atoms needed for one molecule of py- 
roxene, i.e. vg pyr — 1. Since for each molecule we have 
just one Si atom, vg p ■ Asi^pyr = 28. The quantity 
Ai^Pyr = 24 or 56 (magnesium or iron). We must also de- 
fine vfpyr (a number between and 1). Assuming a fixed 
partition x: vf,[ p — 0.8 and Vp^ p ~ 0.2. Finally, 
we have the physical constants for olivines: vg o; = 1, 
Asi,oi = 28, AMg.oi = 24, Ap^^oi = 56, VMg.oi = ^-^^ fi- 
nally Vp^ Qi — 0.4. Since Fpyj. and Fqi are the fractions of 



Mg 



have the following partition of magnesium in dust a 

'^Mg,Ol + '^Mg.Pyr = {^ ■ Fpyr) / [x ■ Fpyr + 2x ■ Foi)a^ig + 

{2x-Foi)/{x-Fpyr + 2x-Foi)a^jg - Algpyra^j^+Mgoia^jg 
and of iron in silicates (once subtracted the iron embedded 
in iron grains) a^^g^, = a^^-a^^p^ = (JFe,oi+'^Fe,Pyr = 

(I-X)- Fpyr/{{1 - X) ■ Fpyr + 2(1 - X) ■ Foi)(T^,^sU + 

2(1 - x) ■ Foi/{{l - x) ■ Fpyr + 2(1 - x) ■ Foi)a^,^sa = 

Fepyrcr^^^szi + F^oi(^Fe,sa- 

It is worth also noting that the dust compounds of which 
we follow the accretion, i.e. pyroxenes/olivines for silicates, 
carbonaceou s grains and i ron grains (the same considered 
bv Izhukovsk a et all (J2008r) ). are fewer in number than the 
dust molecules considered by the stellar yields of dust in 
the most detailed m odels calculated by CNT theory (see 
iPiovan et al.ll2011al for more details). The number and va- 
riety of the dust molecules injected into the ISM are wider 
than the number of accreting species. For this reason, before 
splitting a'^g , dp^ g^i and ag (see below) between pyrox- 
enes/olivines, we need to keep memory of the evolution of 
the other dust species injected by AGB/SNs (also eventu- 
ally swept up and destroyed by shocks in the ISM), which 
are not included in the accretion scheme. 

The Si case. In alternative, silicon can be the key- 
element for olivines and pyroxenes: 



d< 
dt 



— ^ \f Si.Pyr^Si,Pyr ^Si,Pyr] 
+ ^ \f Si,Ol^Si,Ol ~ (^Si.Ol) 



(53) 



where again fsi^pyr and fsi,oi ^^^ ^^^ average condensa- 
tion fractions of sihcon in pyroxenes and olivines. The cor- 
responding partition of silicon in dust is ag — ag p + 

'^Si 01 ~ Fpyr'^si + ^oi'^si- We need to calculate the accre- 
tion time scales of pyroxenes and olivines, to be compared 
with the lifetime of molecular clouds and included into 
the integrations for J Mg,Pyr-, lMg,oh 1 St,Pyr and 1 si,oi- 
Applying Eqn. ((28|) we get: 
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(54) 



(55) 



where 



gr 
^0,Pyr,i 



and 



gr 



are two constants that depend on 



the key-element. Tq p ^ is equal to 2.05, 1.30 and 1.16 for 
the key-elements Si, Mg and Fe respectively, while Tq q; ^ is 
equal to 2.05, 2.60 and 2.32 again for Si, Mg and Fe. A^^ 
and A^p are the atomic weight of the growing species. 
Finally, Ad^oi = Ad^pyr = 121.41. It is important to note 
that in the calculation of the abundances of species avail- 
able for dust accretion, we do not subtract the amount of 
dust already formed a-f from the total ISM abundance a{^ ; 
this is equivalent to assume that the grains of silicates al- 
ready formed act as reactive species in the accretion pro- 
cess. Dropping this assumption (and simplification), the 
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accretion process would be too slow. There is one excep- 
tion though in the case of iron as the key-element for the 
mixture of silicates mixture. In such a case we subtract 
from a{^ the amount of iron already embedded in the iron 
dust grains (that do not react with silicates) , thus obtaining 



_7W _ M 

"i — "Fe 



'Fe,Fe- 



The initial condensation fractions fo,i,Pyr, fo,i,Oi needed 
to derive fi^pyr, fi^oi '^^^ ^^ obtained as in Eqn. ([^ . 

5.3. Oxygen 

Oxygen is much more abundant than the refractory ele- 
ments Si, Fe, Mg, Ca and S. Therefore, oxygen (at least 
the most abundant isotope ^^O) will never become a key- 
element determining the growth of the typical dust grains 
in the ISM. The equation describing the evolution of the 
mass abundance of the ^^O is obtained from Eqn. (|46| . 
dropping the terms giving the contribution of the O as key- 
element. Therefore, considering the accretion of pyroxenes 
and olivines, the temporal evolution of the ^^O abundance 
(we drop the mass number of the isotope) is: 



daS 
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(56) 



where once again the key-element could be Mg, Si or Fe, 
depending on which has the lowest abundance. The_ sub- 
script i indicates a generic key-element. As usual, fi^pyr 

and /j Qi are the average condensation fractions when 
the MC is dispersed, and (jf^py^ and (J^qi ^"^^ defined as 
in Sect. 15.21 Introducing the Dirac delta function nota- 
tion we have: Ai^pyr = 2%5 {A^ - 28) -I- 24<5 [At - 24) -I- 
56(5 {Ai — 56) where Ai is the mass number of the key- 
element i. Furthermore, Ao,Pyr = 16 and, since in pyrox- 
enes we have sulfite SiOs with 3 oxygen atoms, Vq p = 3, 
v[>py^ = 16 [Ai - 28) + 0.86 {Ai - 24) + 0.26 [Ai - 56), and 
■^O.oi = 16. Similarly, Vq qi —4 because in the olivines we 
have sulfate Si04, Ai,oi = 28(5 {A, - 28) + 24<5 (Ai - 24) + 
56(5 {A, - 56), and v^^i = 16 (A, - 28) + 1.6(5 {Ai - 24) -I- 
0.4(5 {Ai — 56). The accretion time scales Tp and Tq^ that 

we need for fiPy^ and fioh ^^e elemental abundances of 
the key-elements and the initial condensation fractions are 
the same as in Sect. 15.21 

In this context, we need also a simple description for 
^^O and ^*0. First of all let us examine how they behave 
in the inert CO molecule. The way in which the isotopes of 
O and C combine to form many isotopologues other than 
12(-i16q jg ygj.y complicated and thoroughly studied, be- 
cause of the importance of CO molecule, its easy detection 
and chemical stability ()van D ishocck & Black 1988; Liszt 
l2007t IVisser et aLll2009D . The photo-dissociation of CO is 
a line process and consequently subject to self-shielding 
that in turn depends on the column density. Therefore, 
isotopologues other than C^^O are not self-shielded un- 
less located very deeply into the molecular clouds (iue to 
the much lower abundance of ^^O and ^^O ()Clavtonll2003t 
iLee et al.l [20081: IVisser et al.]l2009D . Looking at the case of 



the Sun, where the isot( 3pic ratios are about isq/isq g^ 5qq 
and i^o/i^o « 2600 (|Lodders et al.ll20n9t lAsolund et al.l 
2009), the regions with abundances of the isotopologues 
different than C^^O significantly reduced with respect to 
other isotopologues are located very deeply in the atmo- 
sphere. Basing on these considerations, we set the reduction 
factor of ^^O and ^®0, S!co^ equal to 1/3 of the reduction 
factor assumed for CO in Eqn. (P5)) . All this agrees with 
the observational ratios Ci^O/C^O ar id C^^O/C^^O in C 
Oph > 5900 and w 1550 respectively (jSavaee fc Semba(±l 
[l99i). 

Another point to note is that the ratios -'^^O to "'^''O 
and ^^0 in the interstell ar dust could be different from th e 
same ratios in the ISM (|Clavtonlll98llLeshin et al.lll997[ ). 
The subject is a matter of debate, because the above ratios 
strongly depend on th e site where those isotopes condense 
into dust (|Meverj l2009"). Consequently, we should not simply 
scale to ^^O and -'^^O the results obtained for ^^O. However, 
a detailed analysis of this issue is far beyond the purposes 
of this paper and we leave it to future improvements of our 
model. For the time being, basing on the following consid- 
erations we adopt a simple recipe. As in the dust grains, 
"'^^O and ^^0 are much less abundant than other elements, 
consider them as the key element of the dust accretion pro- 
cess. Therefore, Eqn. (|46l) applied to ^^O and -'^^O (briefly 
indicated as "O) becomes: 



da 



— K <^f(^o),Pyr (1 ~ ^Co) '^{"0},Pyr~ 

} (57) 



^0),0l J 



where ^'qq takes into account that not all O and °0 in 
the ISM can be used to form new grains because part of 
them are locked in the CO. 

The relative amounts of ^^O and ^*0 available to form 
silicates are e 



i-o),su = {l-^'co){^Po))/i\-o)0 



the accretion timescales for ^^O and ^^O are: 



T^lTQ 



0.97 



1.00 



"H 
"H 

(JlSQ 



Myr 



Myr 



and 

(58) 
(59) 



where we have considered a typical silicate in which just 
one atom of ^^O and/or ^^O is involved in the accretion 
process. 

5.4. Magnesium 

In our simple picture of the the grain accretion process, 
magnesium intervenes in pyroxenes and olivines, often as 
the key-element. Furthermore, magnesium is present also 
in other dust compounds such as MgO contained in the 
SNffi yields of dust. In analogy to the case of Si, the set 
of equations governing the temporal evolution of the Mg 
abundance is different depending on whether or not Mg is 
the key, similarly to silicon. If Mg is a key element: 
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da 



Mg 



dt 



= K 



{7 



AIg,Pyr'^Mg,Pyr f^M 



K ■ 



(-F M D 

yl Mg,Ol'^Mg,Ol - (^ M g 



Mg,Pyrj + 

,01} (60) 



with the usual meaning of the symbols. The accretion time 
scales and the initial condensation fractions of Mg can be 
obtained in the same way as in Sect. 15.21 following Eqns. 
(jM)) . ([55t and ^^. If Mg is not the key-element, but Si or 
Fe are playing the role, we have the following equation for 
the Mg accretion: 

[ dt )_-''■ A^Pyr^-,,. ^^'.--^»^-- 



D 



n -, , j^ ^Mg,Ol'^Mg,Ol 
<^t,Pyr] + t^ A r.,uD 

■ \J i,Ol"i,Ol ~ "i,OlS 



(61) 



in which AMg,Pyr — 24, ly^/ 

,D 



(sihcon/iron), v^^^py^ 



Alg^Pyr 

1, jy^ 



V 



D 



0.8, A.^pyr = 28 or 56 
= 0.2, AMgM = 48, 
Mg,oi - 1-6, ^S^,Ol = 28, v^^qi = 1 and i/J?^ o; = 0.4. 
The time scales for accretion of pyroxenes/olivines and the 
initial condensation fractions are again from Eqns. (j54p . 
([55)) and dil]). For crf^py^ and af^^i see Sect. 



Fe,Pyr 

28, '^si,oi = 



expression holds good for Aioi- We have also: Ape,Pyr — 
AFe,oi = 56, i^pe.Pyr = 0-2 (according to Mg^Fei-xSiOs 
with a; = 0.8), i^J^py^ = 1-S {A, - 28) + 0.8-6 {A, - 24) + 0.2- 
6 {A, -56), vp^^oi = 0-4 (according to [Mg^Fei-^Ja Si04 
with X = 0.8) and finally v^qi = 1 • 5 {A.^oi - 28) + 1.6 • 
5{A^^oi - 24) + 0.4 • 5{A.,^oi - 56). Also, with the usual 
meaning, f^pyr, f^^oi ^"^^ f Fe,Fe ^-^^ the average conden- 
sations before cloud dispersion. 

For internal consistency, the system of equations must 
take into account that the iron abundance used to describe 
the evolution of the iron grains should be corrected for the 
iron already condensed into pyroxenes and olivines, sub- 
tracting it from apl and (jp^. If the key-element is magne- 
sium, the amount of iron ap^, g^i embedded into pyroxenes 
and olivines is given by: 



D _ f^MgAl 
(^Fe,Sil — ^^ 



^Fe.,Ol , r , ^Fe,Pyr , , 

Mgoi + o ^^9Py' 



D 
Mg,Ol 



V 



Mg.Pyr 



(63) 



where Ap^^ is the atomic weight of the iron and Mgoi and 
Mgpyr are defined in Sect. 15.21 Doing the correct substitu- 
tion we get 

<^Fe.s^l = <^Mg (0-583 ■ Mgoi + 0.583 • Mgpyr) (64) 

In a similar way, if the role of key element is played by 
silicon we have: 



5.5. Iron 

Even considering the small number of accreting compounds 
included in our model, iron can be locked up in grains of 
various type thanks to accretion processes in cold regions 
of the ISM. First of all, iron is locked up in iron grains 
that act as the key element. Second, iron is also present in 
pyroxene and olivine grains. In such a case, iron may or may 
not be the key element. In any case it is removed from the 
gaseous phase and stored in grains. The existence of two 
channels for locking up iron into grains leads to a rather 
complicated equation for the evolution of the iron. There 
are the main terms to consider: in the first and second, 
iron participates to the formation of silicates (pyroxenes 
and olivines) as or not a key-element, and in the third one 
iron plays the role of a key- element for the formation of iron 
grains. Furthermore, the equation will be slightly different 
depending on whether magnesium, silicon or iron is the key- 
element for the formation of pyroxenes/olivines. The final 
equation is 



dt 



r^ ^Fe,PyrVpe,Pyr r-. m 

^ '■ :=: \J i,Pyr"i,Pyr~ 



^i,Pyr'^i^Pyr 



,D 



-CT. 



D X , r^ ^Fe,OlVFe,Ol r-. m 



i,Pyr J 



^^■.oii^F,oi 



or 



~ ^i,Ol) + ^ \f Fe,Fe Y^Fe ~ ^Fe,Sil) ~ 
-{^¥e-<yFe,Sa)} (62) 

where the subscript i stands for Si, Alg or Fe itself. 
Using the Dirac delta function notation, we have Ai^pyr = 
28 -SiA,- 28) + 24:-6{Ai- 24) + 56-S{A,- 56). The same 



'^Fe^S^l = '^S^ (0-8 ' Foi + 0.4 • Fpyr) (65) 



where Fqi and Fpy,. have been defined in Sect. 15.21 If iron 
is the key element for silicates, we simply have: 



'^Fe,Sil — "Fe "Fe.Fe 



(66) 



This clearly requires to keep track of the accre- 
tion/destruction and injection of the iron grains by SNae 
and AGB stars. The accretion time scales for pyrox- 
enes/olivines are the same as in Eqns. (l54)) and ((55|) . 
whereas for iron grains: 



<7H 



'-Fe 



3.5-10-5\ 31.56 



^Fe,Fe 



UH 



{< - - 



D 

Fe,Sil 



Myr. (67) 



where i^pf, p^ — 1 and pp^ = 7.86 gem ■^. Concerning 
the accreting elements we consider the simplest case in 
which they are in form of atoms. The initial condensation 
fractions for magnesium/silicon as key-elements in pyrox- 
enes/olivines follow Eqns. (??) and (??), while for iron we 
have: 



/o. 



Fe,Fe 



X¥e 



.D „D 



'Fe,Sil 



^Fe.max 



M _„D 



0-p,' - cr 



(68) 



Fe.Sil 



where we subtract the amount of iron already locked up in 
silicates and therefore not available to for iron grains. 
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5.6. Calcium 

This element has a comphcate behaviour difficult to fol- 
low. Fir st of all it is usually heavily depleted in the ISM 
(jWhittet 2003; Tielens 2005), and its abundance in the solar 
system is quite low. The difficult arises from the total con- 
densation efficiency and the big fluctuations generated by 
the low abundance, compared to other refractory elements. 
Second there is not an average molecule that could be used 
to represent typical calcium grains in a simple theoretical 
description. Furthermore, the measurements of depletion 
cannot be easily derived from observational data: in prin- 
ciple the ionization equilibrium equation should be solved 
to derive the Ca abundance from the observations of Ca II. 
Possible estimates of Ca abundances by electron densities 
and strengths of the ionizing radiation fields are not easy 
and relying on ratios of ionization and recombination rates 
between different elements (like Cal l to those of Nal or KI 
for in stance) is a cumbersome affair (Weinga rtner fc Draind 
120011) . For all these reasons. Jenkins (2009i, and 2011, pri- 
vate communication) leaved Ca (and also Na and K) aside. 
Trying to overcome this difficulty, we simplify the problem 
as follows. Thanks to its low abundance, we consider Ca as 
the key-element of the associated grains of dust; the equa- 
tion for the evolution of Ca is: 

(^) =KUca,Ca<^^a{Ca)-aSACa)} (69) 

where all the symbols have their usual meaning. To derive 
the accretion time scale t^^ we do not refer to a typical Ca 
dust grain but simply take the shortest timescale among 
those of refractory elements that are most depleted, i.e. 
Mg, Si and Fe {T§l,Tf^ ^Tp^^), and a timescale Tq^^j^ for 
the formation of dust silicates-like molecules with one cal- 
cium atom as key atom. This last-named timescale is mul- 
tiplied for a correction factor Cax, eventually allowing for 
a fast accretion. Therefore t^^ — min{r|[,T^r ,T|.g, Ca;is: • 
'''ca Ca)- '^^^ initial condensation fraction at the MC for- 
mation is as usual jca.Ca = ^cJ'^^a- 

5.7. Sulfur 

Sulfur is a very important element: it is often used as a 
reference case of ne arly zero deple tion in studies of local 
and distant objects (| Jenkins! l2009fl . However the real de- 
pletion efficiency of this element is a matter of debate and 
the assumpt i on of nearly zero depletion cannot be safe. 
ICalura et al.l (|2009[ ) reviews data and theoretical interpre- 
tations gathered over the past years to convincin gly show 
that S can be depleted in considerable amounts. Jenkins 
(|2009f ) points out the depletion of s ulfur can be sig nificant 
along some lines of sight. We take iJenkinsI (|2009[ ) results 
into account to set upper and lower limits to the sulfur 
depletion. Considering sulfur as a refractory element, its 
evolution in dust is given by 



Si and Fe, and a timescale t^s for the formation of dust 



(dal 
\ dt 



) -K{fs^s^P{S)~a^{S)] (70) 

/ accr 

where all the symbols have their usual meaning. Current 
observations do allow us to choose a grain as representative 
of the accretion process. In analogy to what made for cal- 
cium, the accretion time Ts.gr is supposed to be the longest 
between the timescales of the refractory elements, i.e. Mg, 



•s,s 



silicates-like molecules with one sulfur atom as key atom, 
allowing therefore for a slow accretion. We also introduce 
a multiplying scaling factor Sx to eventually correct this 



timescale. Therefore Tg" = Sx ' m&xlTl^jrl^ , T|,g,T|'g}. 
The initial condensation fraction for /g 5 is fs,s " "^ '""'^ 



^^l^i 



5.8. Nitrogen 

Nitro gen is known to be poorly de p leted and if not depleted 
at all (|whittetll2003HTieleni2005[ ). lJenkhIi()2009t) suggests 
that depletion is independent from the line of sight deple- 
tion strength factor and that in any case depletion is very 
low, thus confirming the poor ability of nitrogen to con- 
dense into dust grains. As nitrogen is included in our list of 
elements (both ^^N and ^^N), the evolution of both isotopes 
in dust is governed by: 



Ir) -M7 



„ jvO"iv {N) 



D 



(N)} 



(71) 



where for the accretion time we simply take the longest 
between C and O and the accretion timescale with one ni- 
trogen atom as growing species for a nitrogen molecule. 
Carbon and oxygen are the nearest elements, both do not 
show a strong depletion and, finally, similarly to nitrogen 
both indicate (at least along some lines of sight) low values 
of depletion. Taking the longest timescale, we implicitly as- 
sume that nitrogen very slowly accretes onto dust. As usual 

JN,N — CFnI^N ■ 



6. Supernovae: destruction time scales 

Dust grains in the ISM can be destroyed by other physical 
process such as the passage of shock waves by supernovae 
explosions. The destruction time scale of the element i-h in 
dust grains in the ISM because of the shocks by local SNae 
is defined as the ratio between the available amount of that 
element locked up into dust (the surface mass density of the 
element «-th a^ (fjt), at a given radius and evolutionary 
time) and the rate at which grains containing that element 
are destroyed refueling the gaseous phase : 



TSNR,i = CTj / 



da. 



D 



dt 



(72) 



SNR 



as usual, we drop the dependence from r and t. The de- 
struction time scale can be expressed as follows: 

fdaF 



M: 



destr 



RsN 



(73) 



\ dt y g^j^ 

where Mf*^^*^ = Mf'^'^*^ (r, t) is the amount of mass de- 
stroyed by a single SNa event, while Rsn — Rsn {r, t) is 
the global rate of SNae, obtained by adding together the 
rate of Type I and II SNae. Multiplying numerator and de- 
nominator by ai' and combining together Eqns. (|72p and 
([751) it follows: 



TSNR,; 



ofl 



dal 
dt 



1 



SNR 



Rsn M, 



/f destr 



(74) 



The amount of mass destroyed by the interstellar shocks 
can be defined as the amount of mass swept by the SNa 
shock, multiplied by the fraction of dust mixed with the 
ISM (in this way we get the total swept up mass of dust) 
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and finally by a factor of dust destruction e^ that depends 
on the element we are taking into account: 



M 



destr 



j^ ' ^^^ swept ' ^i • 



Inserting this expression into Eqn. (j74l) we have: 



,-D 



TSNR,'i 



RsNMf^^*^ Msn,eptU\Rs 



T-M 



3N 



(75) 



(76) 



We need now to get an estimate of the mass swept up by a 
SNa shock propagating through the IS M. To describe the 
evol ution of the S Na remnant we follow ICiofR et all i 



and lGibsonI (|1994l ). In brief, the evolution of the remnant af- 
ter the SNa explosion is subdivided into three main phases: 
(i) free expansion that lasts until the mass of the swept 
up material is comparable to the mass of the expelled ma- 
terial; (ii) adiabatic expansion, or so-called Sedov- Taylor 
phase, lasting until when the radiative cooling time of the 
shocked gas is about equal to the expansion time of the rem- 
nant; (iii) radiative expansion, with the formation of a cold 
and dense shell behind the shock front, starting when at 
least some sections of the shocked gas have radiat ed most 
of their thermal energy (lOstriker fc McKegll988t) . Phase 
(i) has a simple solution that can be obtained using the 
obvious relation (4/3) nR^po — Msn where Msn is the 
mass of the expelled material. The duration of the phase 
is r = R/v K, 200 years (very short indeed). For the phase 
(ii) we have t he classical auto-similar a diabatic solution of 
Sedov-Taylor (jOstriker fc McKeelll988ll : 



Rs (i) = 1.15 



En 



Pa it) 



1/5 



t 



2/5 



(77) 



tsf =3.61-10*-£ 



,4 . 3/14,-4/7 







•-H 



where Eq is the energy of the blast wave in units of 10^" erg 
and pg (t) is the density of the environment. The radiative 
cooling leads to the formation of a thin and dense shell at 

^X -5/14 

riH is the number density of hydrogen atoms, Z the metal- 
licity of the ISM and Zq the solar value. Eq is the energy in 
units of lO'^^ erg. The blast wave decelerates until when the 
radiative losses start dominating. At about tpos ~ 0.37ts/, 
the remnant enters the so-called Pressure Driven Snowplow 
phase (PDS) lasting or Phase (iii). In the early stages of the 
PDS phase, the radius of the remnant changes according to 
the following equation (where tsN — t — t'): 

St PDS 3 



Rsit 



SN) 



RpDS 



(79) 



where the radii are in parsec. RpDS is the radius at the 
start of the PDS phase: 



RpDS = 14 • el^'nl^' [ — 



(80) 



The inner side of the shock continuously looses energy be- 
cause of the radiative cooling and it pushes forward the 
shell in the ISM. At the time tmerge, given by: 



5/49 10/49 



Emerge — ZL.l • Tq^Sq Tlf 



z 



15/49 



n) 



The remnant and the ISM loose their identity as single enti- 
ties and merge together. Subsequently, in the time interval 



for tr, 



^ tsN ^ tcooi the radius of the shell is given by: 



R. 



o -D 3/98 3/49 I - 

merge — O.IKpDS^O '^0 I ^~ 

Z0 



9/98 



in parsec, where the time scale tcooi is: 



tr 



203i,f -— 



isf 



Zr. 



-9/14 



(82) 



(83) 



in years. There are four time intervals, with the correspond- 
ing radii, to consider: 



1-15 (;;^j "^ t|jv < tsAT < tpDs; 

RpDS (§71^ - l)"* iPPS < tsN < l.ntsf, 



Rs (tsN) 
Rs (tsN) 

Rs{tSN) =i?PDs(|7^-|)^ 1-1 ''W <iSN <tmerae; 

3 3/ \ -^ 

Rs = i.TRpDSefnf [j^j '" i^erge < tSN < tcooi ■ 

(84) 

To obtain the shock velocity in the various time inter- 
vals, we need the Rs as functions of time. Approximating 
Pg (t) « Pg, we get: 



Vs{tSN) = l^ 
Vs{tSN)^lf^ 



3 tpjjs 
«S (tsN) = 



< tsN < tpDS'y 

tpDS < tsN < 1.17ts/; 

1.17is/ < tsN < tmerge', 
t^merge Ji t'SN ^ t^cool • 



After t 



merge 



(85) 

there is no longer expansion and we simply 
have vs (tsN) = 0. How do we proceed in the chemical 
simulation? First of all we need nn, the number density 
of hydrogen atoms. This is a spatial density, in atoms per 
cubic centimeter; we must derive it from the surface mass 
density of hydrogen in our flat geometry. This is done as- 
suming as parameters the scale height and thickness of 
the disk suggested by the observational data. Once nn is 
(78) known, we derive the time scales tgf, tpjjs and tmerge, and 
finally the radii of the remnant and velocity of the shock 
at the beginning and end of the three evolutionary phases 
we have described. At this point we must fix the velocity 
viow setting the limit at which the shock stops disrupt- 
ing the dust grains . A ty pical value is about 50 km/s (see 
iNozawa et al.ll2006l I2007L for a thorough discussion of this 
point). Knowing the limit velocity, we can soon determine 
the time in which the remnant expanding at the limit ve- 
locity covers a distance equal to its own radius. The swept 
up mass is: 

4 



M. 



swept 



--nR{viow) nn 



(86) 



The other factors of Eqn. (|75)) are the ratio <J^ /Rsn 
and the coefficients e^ that describe the fraction of the el- 
ement z-th condensed into dust that is destroyed by the 
shock. While the ratio a^ /Rsn is simply provided by the 
chemical model in use, to calculate e^ we refer to the studies 
by [No zawa et al. (2006, 2007). In brief, the Noza wa et all 
(|2006t ) data are used in the first evolutionary phases of the 
galaxies to derive the d estruction coef fi cients for the ISM 
dust. From Table (6) of lNozawaet all (|2006f ) we take the 
Qi and bi to be used in: 

e, = a,E'^, (87) 
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where E51 is the energy of the supernova in units of 10^^ 
erg. For simplicity we assume E^i — 1, typical energy of the 
explosion of SNae of mass not too high, that according to 
the IMF are also the most numerous. In this way we get the 
coefficients to be used for destruction by forward sh o cks in 
the ISM. Furthermore, the results bv lNozawaet all ()20Q7[ ) 
are used to constrain the yields by the SNse progenitors, 
before dust is injected into the ISM where it will inter- 
acts with the shocks generated by the surrounding SNae. 
Therefore we take in to account the amou nt of material in- 
jected into the ISM (jNozawa et al.l 120031 ) corrected for the 
effect of the internal reverse shock. 



7. Solar System Abundances 

To validate the results of our model we need a reference 
set of chemical abundances to compare. They are also 
known as the cosmic reference abundances. Usually, the 
Solar System abundances provide the comparison set of 
values, but other choices a re possible (ISavage fc SembachI 
[19961 lAspIund et al.l [20091 lWhittetll2010D . Since the solar 
system abundances can not be determined directly from 
measurements of the ISM, because of the depletion of the 
elements accreted into dust, two main sources of informa- 
tion are usually adopted: spectroscopic data inferred from 
the Sun and Meteoritic Abundances. 

Photospheric abundances of the Sun can not be deter- 
mined directly, but only by means of a model atmosphere 
of the solar spectrum. The model must be able to calculate 
lines formation (classical ID or hydrodynamical 3D), and 
take s into account radiati ve transfer and non-LTE processes 
(see iBasu fc Antial l2008l . for a review about solar abun- 
dances obtained from spectroscopy). However, since there 
are no lines of noble gases in the solar photosphere, coronal 
lines, energetic particles or solar wind are also studied, in 
particular to determine the crucial abundances of He and 
Ne. 

Between the tens of thousands meteorites fell on 
the earth, only fiv e are CI carbonaceous chondrites 
(JLodders et al.l [20091) . They form the only group of mete- 
orites allowing us to determine the abundances in the proto- 
solar environment that match the solar ones, and keep 
some memory of the original proto-environment. However, 
volatile elements and noble gases like H, He, C, N, O, and 
Ne, therefore including the most abundant elements, are 
heavily depleted in these meteorites. We can not rely on 
mass spectrometry to safely determine their abundances, 
but an additional source of information must be adopted. 

It is common practice to combine photospheric 
and meteoritic measurements to get the compilations 
of abundances commonly used in literature and con- 



tinuo u sly updated year after year (Anders & Gr evessd 
1989; Grcv esse fc Noels! 1993; Grcve sse fc Sau val 1998; 
Loddera .20031 lAsplundet al.. .2005a: .Grevesse et al...2007; 



3j 



Lodders et aLil2009HAsplund et al.ll2009[) . 

Two points are worth of mention, (i) According to the 
most recent models, the present photospheric abundances 
of the Sun are lower than the bulk abundances of the proto- 
Sun formed about 4.6 Gyr ago because of the physical 
processes taking place in the deep convective region under 
the solar surface. The effects of thermal diffusion, gravi- 
tational settlir ig due to differential gravity and radiative 
acceleration (jTurcotte fc Wimmer-Schweingruben |2002[ ) 
allow helium and heavy elements to deposit in the interiors 



Table 1. Chemical abundances for the Sun at the present 
age according to the classical compilation by I Grevesse fc Sauvall 
(1998) (GS98 - photospheric s olar a bundances) and the more 
recent ones by lAsplund et al." (2009) (AGS S09 - photospheric 
solar abundances) and Lodders ct alj l|2009D (LPG09 - photo- 
spheric and meteoritic abundances compilation). Finally, in the 
last column we present the abundances we have adopted in this 
study. The Abundance A(X) of the element X is in units of 
logj^Q (N(X)/N(H))-|-12. Only the elements included in our chem- 
ical model are reported and discussed 



Element 



GS98 



AGSS09 



LPG09 This work 



H 


12.00 


12.00 


12.00 


12.00 


He 


10.93±0.004 


10.93±0.01 


10.93±0.02 


10.93±0.02 


C 


8.52±0.06 


8.43±0.05 


8.39±0.Q4i 


8.50±0.06^ 


N 


7.92±0.06 


7.83±0.05 


7.86±0.122 


7.83±0.05 





8.83±Q.06 


8.69±Q.05 


8.73±0.07=^ 


8.73±0.06 


Ne 


8.08±0.06 


7.93±0.10 


8.05±0.10'' 


7.99±0.10 


Mg 


7.58±0.05 


7.60±0.04 


7.54±0.065 


7.57±0.05 


Si 


7.55±0.05 


7.51±Q.03 


7.53±0.Q1'^ 


7.53±0.01 


s 


7.33±0.11 


7.12±0.03 


7.16±0.02^ 


7.14±0.03 


Ca 


6.36±0.02 


6.34±0.04 


6.31±0.02'' 


6.34±0.04 


Fe 


7.50±0.05 


7.50±0.04 


7.46±0.08^ 


7.50±0.04 



z/x 0.0231±0.018 



^This is the same low photosphere value as in iLodderi (|2003l ), 
selected fromi Allcndc Prieto et al.| (|200^ ), an d confirmed in 
lAsplundetaL, (2005b) and lScott et al.l l|2006l ).^N is taken 
from lCaffau et al. (200£|), from solar photospheric 

models. ^The O abundance is an average from solar 

photo s pheric mode ls by Caffau ct al. (2008); Ludwi g fc SteffenI 
(|2008D : lMeiend ez fc Asplund (2008). "Nc a bundan ce is an 
averag e from [Morel fc Butler (2008); Lan di et all 
(|2007 ). ^Average between solar and me t eoriti c 



values. "Meteoritic value. ICaffau et aD (|201Gl ) 



of the Sun. The deca y of radioactive elements affects the 
isotopic compositions (jPiersanti et al.l [20071 : iLodders et al.l 
120095- The correction needed to obtain the unfractioned 
abundances of the proto-Sun is of about 0.05 dex for 
He and 0.04 dex for heavier elements (jAsplund et al.l 
2009) o r slight l y hig her of 0.01 dex for both corrections 
(Lodders et al.l l2009l) . With these corrections there is 
ag reement with th e abun dances observed in nearby B stars 
bv lPrzvbillaet al.l (|2008[ ). (ii) The He and Ne abundances 
determined with coronal lines or solar wind suffer from 
the FIP effect. Therefore the most precise method to 
determine the He abundance is helioseismology and the 
adiabatic index Fi, whereas for Ne abundance are the 
photospheric ratios_Ne/0 and Nc/Mg together with their 

systematic errors (Ba su & Antia 2008)_; 

Starting fronr JAllen de Priet o et aD (|200lL I2002D . 
the widely used compilation of abundances by 
iGrevesse fc Sauval (JH)98) has been the subject of a 
continuous revision of the photospheric C, N, O abun- 
dances toward a significantly lower value, i.e. abo ut 0.2 dex 
lower than before (| Asplund et Id] l2004l l2005al[b[ ). These 
new values, derived from 3D hydrodynamical models with 
updated input physics and no macro- and micro-turbulence, 
bring the Sun to better a gree with the metal conten t of the 
Galactic neighbourhood (JTurck-Chieze et al.ll2004l) . They 
are also in good agr eement with recent determ inations of 
B-star abundances (jNieva fc Przvbillal l2008bl [al) and car- 
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bonaceous CI chondrites fo r many elements (lAsplund et alJ 
l2005bt iLodders et al.ll2009t lAsplund et al.ll2009D . There is 
also a better agreement between different diagnostics (01, 
[01] and OH vibrational and roto-vibrational bands for 
Oxygen, CI, [CI], CH, C2 and CO for Carbon), but a 
strong disagreement with the helioseism o logical models 
(Antia k Basu"2005'; 'Bah call et alll2005allbl : lAntia fc Basiil 
1.2006 : Basu fc Antia 200j). In brief, the position of the 
base of the convective zone, the helium abundance, sound 
speed and d ensity profiles are a ll affected by the revised 
abundances JBasu fc Antial 120081 ) . Many suggestions have 
been advanced to solve this discrepancy, e.g. increased 
input opacities, increased abundances of neon or other 
elements, fine tuning of the diffusion process, inclusion of 
other additional physical processes etc. All this, however, 
without removing the discrepanc y. The recent r evisio n 
of the abundance compilation by lAsplund et al.l ()2009() . 
reversed the trend toward higher abundances, thus in 
closer agreement with helioseismology, even if the dif- 
ferences are still significant. Also, a series of models 
calculated over the past few years (all based on the 
so-called CO^BOLD 3D hydrodynamic al code) go toward 
a reconciliation with h e Uoseismolo gy (ICa ffau et al.l l2008t 
Caffau fc Ludwid 120081: ICaflFau et al 2009 LlLudwig et al.l 



20091 iMaiorca er al."2009': 'Caffau et al.ll2010li . Anvwav. as 
pointed out in [Caffau et al. (201C|), results from different 
3D hydrodynamical simulations may still mutually differ 
by as much as 0.1 dex, suggesting that better validation of 
the hydrodynamical mo dels and more and more updated 
solar models are needed (jTurck-Chieze et al]|2008[) . 

In Table |T] we sh o w the classical set of abundances 
by iGrevesse fc Sauvall (|1998l ) (GS98) the two most recent 
compilations availabl e in literature b v lAsplund et al.l (|2009[ ) 
(AGSS09) and Lodd ers et ail (l200l (LPG09) and, finally, 
in the last column, our adopted reference set for the present- 
day solar abundances. The abundances by number of the 
various elements are indicated by A(X) which stands for 
logio(N(X)/N(H)) + 12.^ 

Let us now examine in some detail the choice we made 
for the abundances of the various elements included in our 
model: 

Carbo n: the carbon abundance is A(C)=8.50±0.06 ac- 
cording to ICaffau et al.l (|2010[ ). It corresponds to a weak 
efficiency of the collisions with neutral hydrogen atoms 
(Sh — 1/3). The uncertainty is the sum of statistical (0.02) 
and systematic (0.04) errors. This value is similar to the old 
GS98 determination and very close to the average estimate 
of the carb on abun dance during the past thirty years (see 
Table 4 in Caffau et al.ii20ld) . 

Nitrogen: the abundance of Nitrogen determined by 
AGSS09 does not much differ from that by LPG09, but 
it is significantly lower than GS98. We adopt the value 
by AGSS09 A(N)=7.83±0.05. As this is derived from both 
atomic and molecular lines, the error is minimized. 

Oxygen: Oxygen, together wit h Neon, are the mos t 
debated and uncertain elements (lAspl und et al.l l2009fl . 
LPG09 adopted the average value be tween the estimate 
from CO^BOLD at mosphere models bylCaffau et al.l (|2008D 
(8.76±0.07), tha t bvlLudwig fc SteffenI (I2008D (8.72^0.06'!. 
and the one bv iMelendez fc A splund (2008) (8.71±0.02). 
AGSS09 recommends 8.69±0.05 obtained from atmosphere 
models with updated physics and taking the mean value 
of estimates based on [01], 01, O H vibration-rotation and 
OH pure rotation. IWhitted (|2010D adopts 8.76±0.03 based 



upon observations of nearby B-stars (jPrzvbilla et al.ll2008l ). 
Averaging all these values, we adopt 8.73±0.06. 

Neon: the abundance of Neon proposed by AGSS09 
(7.93±0.10) is determined from their O abundance and the 
Ne/0 ratio. LPG09 (8.05±0.10) takes the mean value be- 
tween measurements in B-type stars and UV-flares. The 
proposed values are very different. We adopt the aver- 
age value of 7.99±0.10. The corresp onding Ne/0 ratio 
is 0.182 in ver y goo d agreement with lYound ()2005D and 
ISchmelzetaD ((20051) . 

Magnesium: magnesium ^/-values for the two ioniza- 
tion stages Mgl and Mgll suffer of well know uncertainties, 
as underlined in AGSS09. For this reason we adopt the 
weighed average between the LPG09 meteoritic abundance 
7.55±0.01 and the AGSS09 value, obtaining 7.55±0.05. 

Silicon: silicon is the element linking meteoritic and 
photospheric measurements. Since H is very depleted in CI 
chondrites, meteoritic abundances are usually expressed in 
the cosmochemical scale as number of atoms per 10® Si 
atoms. If the abundance by number of Si is known, we may 
easily obtain the abundance by number of any element in 
the usual scale A(X) = log^o (N(X)/N(H)) + 12.' The various 
estimates of the Si abundance listed in Table [1] agree quite 
well each other. Therefore we can adopt the estimate by 
LPG09 that is the average between the solar photospheric 
and meteoritic values. 

Sulfur: the sulfur abundance has been revised by 
AGSS09, including additional lines and NLTE corrections. 
We take the mean between this value and the meteoritic 
abundance by LPG09 (see their Table 4). 

Calcium: for the abundance of Calcium we choose the 
value proposed by AGSS09, which is in good agreement 
with LPG09 and GS98 a.nd is based upon recent NLTE 
abundance corrections bv iMashonkina et al.l (|2007l ). 

Iron: we adopt the value given by AGSS09, in excellent 
agreement with the meteoritic measurement by LPG09 and 
GS98. 

Once the reference abundances are assigned, we can 
tackle the element depletion in the ISM, that is the under- 
abundance of some elements compared to the correspond- 
ing reference value, because a fraction of the element under 
consideration is locked up in the interstellar dust. The de- 
pletion of the element X in the ISM i s measured in the fol- 
lowing way. Following iJenkinsI ()2009[) and Whittct (201(it ) 
we define 



[Xgas/H] = log [N (X) /N (H)] - log (X/H) 



ISM 



(88) 



where N (X) is the column density of the element X and 
N (H) = N (HI) +N (H2) the same for hydrogen and, finally, 
and (X/H)jgj^ the assumed reference abundance. If part of 
the element X is locked up into dust grains ([Xgas/H] < 0), 
the fractional abundance of the element depleted into dust 



(Xdust/H) = (X/H) 



ISM 



1 - 10'^' 



[X.^JH] 



(89) 



Usually, establishing the degree of depletion of an el- 
ement is a cumbersome affair, because the choice of the 
reference abundance is a difficult task and errors and in- 
consisten cies affec t the m easurements (Whittct 2003, 201(|). 
Recently, [Jenkins! (|2009D presented a thorough study of the 
depletion for 17 elements along 243 lines of sight, trying to 
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focus on the rates of depletion and leaving aside the prob- 
lem of the reference abundance. He was able to characterize 
the degree of depletion according to three parameters: one 
describes the overall level of depletion along a particular 
line of sight; the other two are related to the depletion of 
each element along that line. Once ado pted the reference 
abu ndances proposed by 'Loddcrs' (2003) for the proto-sun, 
iJenkins (2009 ) presents the depletion parameters for the 
smallest and largest depletion efficiencies (his parameter 
F*), thus bracketing the region to be matched by any the- 
oretical model. In the context of the present study, to cor- 
rectly make use of the 'Jenkins' (2009) results, the following 
rem ark is appropri ate. Jenkin s (2009) analysis stems from 
the iLoddersI ()2003D present reference abundances (see his 
Table 1) which is different from the one we have adopted 
(see the entries of Table [1] Theref ore, a su i table shift 
has to be applied. Furthermore, the ILoddersI ()2003D ref- 
erence abundances for the proto-solar environment have 
been obtained from the present-day solar -|-meteor itic ones 
by a dding ~ 0.07, which is likely too a correction (" Jenkins! 
[2009 : Loddcrs ct al. 2009; Przvbilla ct al. 2008). Therefore, 
we derive our present-day reference abund ances applying 
the m ore realistic correction suggested by (jAsplund et al.l 
[2OOI . 

Table [2] shows the reference abundances for the SoNe 
used by Jenkins (200!|) and those adopted in our set for 
the proto-solar environment, together with the upper and 
lower limits for the smallest (F* = ) and heaviest ( F^ = 1) 
depletion efficiencies adapted from iJenkinsI ()2009f ) to our 
case. Calcium deserves some remarks. It is not included 
in the I Jenkins! (|2009il list, whereas we take it into account. 
Furthermore, as this element is heavily depleted in the ISM 
(|Whittet ll2003l: iTielenf^ 12005"), we consider the values pro- 
posed by Crinklaw ct al. ( 1994) for the smallest and h eav- 
iest depletion. The value proposed by !Whittet! ()2003!) for 
the diffuse clouds falls in the middle. 

In Fig. [2] we finally show the range of depletions adopted 
in our work to test the simulations versus the SoNe. Since 
the depletion is line-of-sight dependent, we show for each 
element the range allowed taking into account the smallest 
and largest depletion with relative error bars. The values 
of depletion for cool and warm disk are also shown. For the 
sake of clarity, we expanded the depletion range for C, N, 
O in the lower left corner of the figure. 

8. Models of the ISIVI with dust: results 

The chemical model for disk galaxies we are proposing ex- 
tends the the original dust-free, multi-zone model with ra- 
dial flows de veloped by the Pa dova group over the years 
(!Chiosi"198rt 'Chfosi fc Maeder"1986': 'Portinari et al.!|1998!: 
[POTtinari & Chiosi 1999, 2000; Portinari et al. 20043^ to 
which the reader should refer for all the details not men- 
tioned here. The model is quite complicate and obviously 
contains many parameters. Therefore it would be wise to 
suitably select the parameters to vary guided by some gen- 
eral considerations to be kept in mind: (i) This study is 
mainly devoted to highlight the role of dust in chemi- 
cal models rather than perfectly reproducing the proper- 
ties of gas and stars of the MW in the local pool. Again 
dust properties will be the target of a forthcoming paper, 
just as fun ction of the galacto- centric distance along the 
MW Disk (| Piovan et al.! !2011c! ) . For this reason various 
recipes for dust formation and evolution in the ISM must 
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Fig. 2. The range of depletions is shown for those ele- 
ments in our list that intervene in the dust formation or 
are present in the ejecta by AGB stars and SNa; (C, N, 
O, Ms, Si, S, Ca and Fe). The insert shows in detail the 
depletions for C, N and O. Four values are plotted for each 
element. Filled squares and circles represent the smallest 
and largest depletion, respectively. The crosses determine 
the range for the depletion in the War m and Cool com po- 
nents of the Galactic Disk , according to !Jenkirisl (|2009! ) and 
ISavaee fc SembachI (|1996! ). See Table [2] for more details. 

be tested, (ii) A reasonable agreement between model re- 
sults and observational data for gas, dust and stars (both 
in the solar vicinity and across the galactic disk) must 
be reached for the sake of physical consistency of the 
whole model, (iii) Finally, we will take advantage of the re- 
sults a l ready obtained by Chios i (1980); Chiosi fc Maederi 
(119861 ) : !Portinari et a l. (1998): !Portinari fc Chiosil (119991 
I2OOO! ): !Portinari et al.l (|2004bJ ar) as far as the fine tuning 
of several important parameters is concerned. 

In the following each model is identified by a string of 
nine letters (the number of parameters) in italic face whose 
position in the string and the alphabet corresponds to a 
particular parameter and choice for it. The position in the 
string is the same as in the list below. Let's now shortly 
comment on the parameters we have considered and the 
choices we have made for each of them together with the 
identification code. 

(1) The IMF with its lower and upper limits and the frac- 
tion C of stars with mass M > IMq. Eight IMF are 
considered as described in Sect.|31 Salpeter (A), Larson 
(B), Kennicutt (C), Kroupa original {V), Chabrier {£), 
Arimoto (J"), Kroupa 2002-2007 {G), Scalo {U) and, fi- 
nally, Larson adapted to the SoNe (J). Since some of 
them are similar, we will examine in particular only the 
results obtained for some interesting cases useful to un- 
derstand the influence of the IMF. The Upper and lower 
mass limits (and Q are selected according to the default 
values already given in Sect. [3] 

(2) The Star Formation law. Five SFR are considered (Sect. 
H constant SFR (A), Schmidt law (B), Talbot fc Arnett 

(TfcA - C), Dopita fc Ryder (DfcR - V) and, flnally, 
Wyse fc Silk (WfcS - £). 
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Table 2. In column (1) we list the eleme nts we have con sidered. In columns (2) a nd (3) we show the chemical abundances for the 
proto-solar environment, according to the lLoddersI ()2003l ) compilation adopted by iJenkina ()2009l ) and the values we have adopted. 
In columns (4 ) and (5 ) we g ive [Xgas/HJ^ and [Xgas/H]j^, the smallest and largest depletion values for the various elements are 
obtained from iJenkinsI (|2009l ). with a correction to account for the difference between the reference set of elemental abundances. 
In column (6) we list the abundances in the ISM and in columns (7) and (8) we show the abundances of elements in dust in units 
of 10® (Xdust/H), for the smallest and largest depletions, res pectively . Finally colu mns (9) and (10) sh ow the abundances for the 
Warm Disk and the Cool Disk, according to the fit made by IJenkina l(200S' l of the lSavaee fc SembachI (|199Q) data 



Element A{X)'^ 



A{Xf 



H 

(■10® 



(■10® 



(■10® 



CD 



(1) 



(2) 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



(9) 



(10) 



H 


12.00 


12.00 


- 


- 


10® 


- 


- 


- 


- 


He 


10.984±0.02 


10.98±0.02 


- 


- 


9.5^10* 


- 


- 


- 


- 


C 


8.46±0.04 


8.54±0.06 


-0.192±0.194 


-0.293±0.075 


347 


124t?°, 


1701^3 


-0.204 


-0.283 


N 


7.90±0.11 


7.87±0.05 


-0.079±0.119 


-0.079±0.119 


74 


I2tll 


12li^ 


-0.079 


-0.079 


O 


8.76±0.05 


8.77±0.05 


-0.020±0.060 


-0.246±0.055 


588 


27lg 


255j;« 


-0.047 


-0.222 


Ne 


8.08±0.06 


7.97±0.10 


- 


- 


93 


- 


- 


- 


- 


Mg 


7.62±0.02 


7.61±0.05 


-0.260±0.030 


-1.257±0.029 


41 


18+^ 


38+1 


-0.380 


-1.157 


Si 


7.61±0.02 


7.57±0.01 


-0.180±0.035 


-1.319±0.052 


37 


13+2 


35+1 


-0.320 


-1.205 


S 


7.26±0.04 


7.18±0.03 


0.243±0.092 


-0.635±0.206 


15 


-ntr 


utl 


0.137 


-0.548 


Ca 


6.41±0.03 


6.38±0.04 


-2.090±0.200 


-3.930±0.200 


2.4 


2.38 


2.398 


-2.311 


-3.746 


Fe 


7.54±0.03 


7.54±0.04 


-0.951±0.038 


-2.236±0.041 


35 


3lti 


34t° 


-1.105 


-2.107 



Abun dances A{X) of the element X are in units of logjQ (N(X)/N(H)) + 12. ^Proto-solar abundances adopted in (|Jenkinsl 
120091 ). ^Proto-solar abundances adopted in this work. '^The positive depletion of sulfur means that instead of observing a 
depletion of the element the gas seems to be enriched of sulfur atoms respect to the reference set. 
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Fig. 3. Temporal evolution of the contribution to the abundance of dust by the four types of grain (on which we distributed 
the single elements) and the three sources. All the contributions have been properly corrected for the destruction of dust. 
Left panel: results for an inner ring of the MW centered at 2.3 kpc. We show: silicates (continuous lines), carbonaceous 
grains (dashed lines), iron dust (dotted lines) and, finally, other grains bearing S, Ca and N (dot-dashed lines). For 
each group we distinguish the net contributions from the ISM accretion, AGB and SNse, that is: ISM-C, AGB-C and 
SNffi-C for carbon grains, ISM-Sil, AGB-Sil and SNae-Sil for silicates, ISM-Fe, AGB-Fe and SNa;-Fe for the iron dust and 
finally, ISM-Ca/S/N, AGB-Ca/S/N and SNffi-Ca/S/N for the other grains. The two thin dashed lines represent the same 
GVABBCBBB model but with 15% (upper dashed line) and 45% (lower dashed line) of CO. In all cases ^co = 0.30 as 
default value. Central panel: the same as in the left panel but for the SoNe at 8.5 kpc. Right panel: the same as in 
the left panel but for an outer ring at 15.1 kpc. Since the ratio (Td/o'h is lower, the scale of the y-axis is shifted respect 
to the scales for the inner region and the SoNe in the left and central panels. 
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(3) 



(4) 



(5) 



(6) 



The fraction of MCs in which dust accretion takes place. 
This quantity must be specified because the chemical 
model in use does not contain a real multi-phase de- 
scription of the ISM. Two cases are included: a constant 
fraction (A) based on the SoNe data {xmc — 0.2) and 
a varying xmc (S) related to the local SFR and total 
gas density (see Sect. 14. 2|) . 

The model for the accretion of grains in cold molecu- 
lar regions . Two choices are available: the simple model 
by IDwekl (|1998[ ) based on ty pical accretion times cales 
for dust grains as modified by ICalura et al.l ()2008D (see 
Sect. 14.11 (A) and t he rec ent and more refined model 
by Zhukovska et al.l (|2008[ ) (see Sect. SI]) (S). 

The condensation efficiencies for d ust in Type la 

SNae. Two cases are p ossible: iDwekl ()1998[ ) (A) and 
IZhukovska et al.l (|2008h (B). In the former they con- 
tribute to the dust budget in a away comparable to that 
of Type II SNae, whereas in the latter their role is neg- 
ligible except for a small amount of iron in agreement 
with the observations. The diff erent choices f o r this 
rameter have been discussed in lPiovan et al.l (|201l£ 

The condensation efficiencies for dust in Type I I SNae . 
Three choices are possible. The first one bv lDwekl ()1998[ ) 
who suggest s a high condensation efficiency ( A), the ef- 
ficiencies bv lNozawaet al.l (|2003l I2006L l2007tl based on 
dust nucleation models and taking into account the ef- 
fects of t he reverse and fo r ward shocks (B), and finally 
those bv IZhukovska et aTl (|2008l ) who favor a low con- 



(7) 



(8) 



densation efficiency in SNae based upon pre-solar grain 
observations. The impact of these dif ferent choices o n 
the results has been discussed in , Piovan et al.l (|2011al ). 

The condensation efficiencies for dust from AGB stars. 
Two c hoices are available: the simple recipes by IDwekl 
(|1998| ) (A ) and the condensation efficiencies obtained 
by full cal culations of dust format ion in synthetic AGB 
models bv lFerrarotti fc Gail (|2006D (B). Again, for more 
details on these different pos s ibilities for AGB Stardust 
production, see lPiovan et al.l (|2011aD . 

The age at which a bar is introduced to reproduce the 
radial distribution of the gas in the Galactic disk, in par- 
ticular in the region of the molecular ring around 4 kpc. 
For the purposes of this work we do not play with the 
bar effect: the pattern of velocity to simulate the bar ef- 
fect in the radial flows me chanism is simply taken from 
iPortinari fc Chic)sil ([2000), (to which the reader should 
refer for all the details) and is suitably chosen for ev- 
ery SFR law. There is no bar-effect on the SoNe, for 
this reason between the three cases included (no bar ef- 
fect (A) , an onset of the bar 4 Gyr ago (B) and 1 Gyr 
ago (C), we will simply just fix the bar effect at the 
cas e (B). For more deta ils and discussion on this point, 
see lPiovan et al.l (|201ld ). where the dust formation and 
evolution on the whole disk is examined. 

The parameters of the SFR l aws: these are chosen ac - 
cording to the discussion by (jPortinari fc Chiosilll999D . 
While the exponents fc or m and n are fixed (see Sect. 
|3|) we let the efficiency of the star formation v vary and 
ass ume three value s more or less in the ranges suggested 
by (jPortinari fc Chiosi 1999) for every SF law. The min- 
imum value for Vmin is case A and the maximum value 
for Vmax is case C). The average case is B. 



Table 131 summarizes the parameters we have just described 
together with their associated identification code: the se- 
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Fig. 4. Temporal evolution of the contribution to the abun- 
dance of dust during the first 1.5 Gyr- 2 Gyr. All the con- 
tributions have been properly corrected for the destruction 
of dust. Left panel: results for the inner ring of the MW 
at 2.3 kpc. We show for the SNae injection: silicates (thin 
continuous line SNae-Sil) , carbonaceous grains (thin dashed 
line, SNffi-C), iron dust (thin dotted line, SNae-Fe) and, fi- 
nally, other grains bearing S, Ca and N (thin dot-dashed 
line, SNag-Ca/S/N). The thick lines represent the total con- 
tribution from one source to the dust budget sub-divided 
in: SNae (thick dashed line, SNae-tot), ISM (thick continu- 
ous line) and AGB (thick dotted line). Central panel: the 
same as in the left panel but for the SoNe at 8.5 kpc. Right 
panel: the same as in the left panel but for the outer ring 
at 15.1 kpc. 



quence must be read from top to bottom. For example, the 
string VBAABABAB corresponds to Kroupa 199 8 IMF, 
Schmidt SFR, ANN model for ^_mc, IDwekl (|1998D accre- 



tion model, IZhukovska et al.l (|2Q08f ) type la SNae recipe for 
dusty yi elds, Dwek (11998 ) c onden sation efficiencies for type 
II SNaj, iFerrarotti fc Gaill (|2006D yields for AGB stars, no 
bar and high efficiency v of the SFR. If not otherwise spec- 
ified radial fiows and bar effect will always be included by 
default. 

Obviously, these are not the only parameters of chem- 
ical models: even the classical ones are themselves very 
rich of parameters and when including also radial flows, 
bars and dust formation/evolution the parameter space ac- 
quires many more dimensions so that a full exploration of 
it is a cumbersome affair. However, it may happen that: 
(1) some variables play a secondary role, thus not influenc- 
ing that much the results (in particular we focus on the 
effects related to dust, the target of the work) and just 
adding second order corrections; (2) some variables, even 
if playing an important role and influencing signiflcantly 
the results, have a well deflned and/or restricted sphere of 
influence therefore are not of much interest here because 
their effect is clear. The parameters belonging to these two 
classes are not varied in the models, but kept flxed to a 
suitable value. In the following we examine the effect of 
some primary and secondary parameters and the general 
behaviour of the dust model. To this purpose we adopt the 
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Table 3. Parameters of the models. Column (1) is the parameter number, column (2) the associated physical quantity, and 
column (3) the sourc e and the italic symbols are the identification code we have adopted. See the text for some more details and 
iPiovan et al.l (|2011bt ) for a detailed description. 



Parameter 



Source and identification label 



IMF 

SFR law 

Xmc model 

Acer, model 

SNae la model 

SNae II model 



7 AGB model 

8 Galactic Bar^'' 

9 Efficiency SFR^'^ 



Salpeter^ {A), Larson^ (B), Kennicutt^ (C) Kroupa orig.* ("D), 

Chabrier^ {£), Arimoto** {T), Kroupa 2007'' {G), Scalo* {W), Larson SN^ (X) 

Constant SFR {A), Schmidt^^ (S), Talbot & Arnett" (C), Dopita & Ryder^^ (D), Wyse & Silk" 

Artificial Neura l Netw orks model^^ (^), Const ant xmc as in the Solar Neigh. ^ ^ {B) 

Modified IDw3 (lf998l) andlCalura et al.' (2008') {A); adapted 'Zhukovska et al." f200^) model (B) 

Dust injection adapted from: Pwck (1998), Calura et al. (2008) {A), Zhukovska ct all l|2008l ') (B) 

Dust injection adapted from: iDwekl (|f998l ) f.4). IZhukovska et all (|200al (S), 

iNozawa et all (|2003l . [2006I . [2007^ (C ) 

Dust injection adapted from: iDwekl (|1998l ) (.4) . iFerrarotti fc Gaill (|2006H {B) 
No onset [A) , onset at tc — 4 Gyr (B) , onset at tc — f Gyr (C) 
Low efficiency (A), medium efficiency (B) , high efhcicncy (C) 



{£) 



Saloeted (Il955l'l. ^Larsonl lll986l 11998 



(12001 1. '^Arimoto fcYoshiil ||1987D. iKroupal ('2002b, 2007J) 

(2004b). "'Schmi dtlTigSgfj. ^HTalbot fc Arn ctt (1975J. "Tal bot fc Arnet .__ _. . ____ ,_ 

^ jPiovan et al.i (|2011d ). ' iZhukovska et alf (|2008i V ^ IPortinari fc Chiosil (120001 ). ^ iPiovan et 311(201 Ibl 



iKennicutt et al.l (119941 ). iKroupal (119981). iChabriei 

(Il986f). Uarsonl (119861) : iScalol (119861): iPortinari et af 

'19751) . 'JDopita fc R vdcr ( 1991). ' ^Wvse fc Silkl (|198 ^ 



QVABBCBBB m odel as the 'default' one, with radial flows 
and bar included (jPortinari fc Chiosill2000t) as the reference 
case. This model uses the most detailed theoretical recipes 
for the yields of dust and the amounts of MCs, while the 
SFR and the IMF are simply selected between the available 



ones. 



8.1. Fraction of CO in t/)e ISM 

As already discussed in Sects. [01 and l5.3l a percentage from 
20% to 40% of the Carbon and u p to 20% of the Oxygen is 
locked in the CO molecules (see Ivan Dishoeck et aLlll993t 
Ivan Dishoeck fc Blakd Il998l ) . This affects the amount of 
Carbon available as key-element for the growth process ow- 
ing to the lowered abundance of this element in the gaseous 
phase. Furthermore, in the dense and cold reg ions of the 
ISM CO tends to condense onto d ust gr ains (jGoldsmithl 
120011: iBacmann et~alll2002t IWhittetl [20Tol). leaving th e gas 
phase and sticking efficiently to dust (|Whittetll2010t) . The 
depletion of the CO can reach in pre-stellar cores ev en a 
factor of 10 (|Tafalla et al.l[200llWah^slev et al.l[200l . All 
this is parameterized by the fraction C^co of CO. The effect 
of Ceo on the key-elements ^^C and ^^C (and of Ceo o'^ ^^^ 
and ^^O) is straightforward: the higher is the amount in 
CO, the lower is the amount of free atoms in the gas-phase 
available for the growth, thus implying longer timescales 
and slower dust formation. Therefore, we expect a lower 
contribution by the accretion to the dust budget of the 
ISM when high values of Ceo are adopted. 

In Fig. 121 we present the time evolution of the contri- 
butions to the dust budget for three regions of the MW: 
a central region (r^ = 2.3kpc), the SoNe (rfc = 8.5kpc) 
and an outer region (r^ = 15.1kpc). First of all, we do 
not show the temporal evolution of the single elements, 
but of important groups of elements representative of 
the main typical ISM dust types. These are the silicates 
(olivines-|-pyroxenes-(-quartz-|-silicon in SiC), carbonaceous 
grains (carbon grains-|-carbon in SiC), iron grains and, fi- 
nally, other grains containing S, N and Ca. All the dust 
abundances have been normalized to the local hydrogen 



density and corrected for dust destruction so as to repre- 
sent the effective net contribution to the dust budget. The 
contributions are split in three main sources: AGBs, SNae 
and ISM accretion. In order to explore the effect of CO, in 
the QVABBCBBB model we vary the CO abundance from 
15% up to 45%: the two thin dashed lines in Fig. |3| bracket 
the region of variation of the ISM contribution. Carbon ac- 
cretion in the ISM is the only process varying with the CO 
abundance: the lower ^co the higher is the ISM contribu- 
tion. This effect can be significant and we cannot easily 
get rid of it when evaluating carbon depletion: however, it 
is clear and straightforward and limited to the only car- 
bon. Even if ^^O is a component of the CO molecule, due 
to its high abundance and never being a key-element, it is 
scarcely affected by variations of Ceo • The same holds true 
for ^^O and ^*0 but in this case due to their low abun- 
dance: their budget depends on ^co , but their contribution 
to the global budget is negligible. Therefore, the variations 
of Ceo do not affect the silicates budget. 

From the analysis of Fig. |3| we can notice several gen- 
eral features of the dust evolution common to all the models 
that are worth to be underlined: (i) the main contribution 
to the dust enrichment during most of the Galaxy lifetime 
is due to the accretion process in the ISM. Furthermore, 
dust production is much higher in the inner regions of the 
Galactic Disk compared to the outer ones where the weight 
of the ISM gets smaller so that during many Gyr the star- 
dust injected from AGB stars and SNae drives the total 
yields. This is ultimately due to the low number densities 
of metals that do not favour the accretion process; (ii) the 
inner regions reach higher metallicities than the outer ones: 
therefore, all physical processes depending on the metallici- 
ties are much enhanced in the central regions of the Galaxy, 
e. g. the yields of dust from oxygen rich M-stars; (iii) even 
if some sources (like AGB-C stars in the inner regions) 
may vary their contributions with time, in general the total 
amount of dust keeps growing monotonically. In some way, 
this mirrors the metallicity enrichment of the ISM and the 
fact that dust formation is very sensitive to metals. 

In Fig. |4|we show in detail the evolution of the contri- 
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butions of AGB stars, SNa^ and ISM to the total budget of 
dust during the first 1.5-2 Gyr for the three selected regions 
of the Galactic disk. For AGB, SNae and ISM we show the 
total budget. However, limited to SNs, we also distinguish 
the various types of grains. During the early stages, star- 
dust dominates the scene: going from the innermost to the 
outermost regions it takes more and more time for the ISM 
accretion to overcome the stellar contribution. This is ul- 
timately due to the higher densities and metal content of 
the inner regions that favors the onset of accretion. For the 
outer regions, with low density and low SFR, the ISM starts 
to be important only at t>4 Gyr (see Fig. [3]). It is interest- 
ing to note how the dust-to-gas ratio for net yield by SNae is 
nearly constant going from the innermost regions (left pan- 
els of Fig. U and 131) to the external ones (right panels of the 
same figures), even if the associated timescales are much 
different. Toward the center of the Galaxy, we have more 
hydrogen, higher SFR, higher SNk rate and higher yields 
of dust. Both gas (via infall) and SNae dust (via infall, SFR 
and SNae explosions) grow in absolute value moving from 
outside to inside, whereas their ratio does not change in the 
same proportions (considered that in any case the yields of 
dust by SNae do not depend on the metallicity, at least ac- 
cording the the kind of theoretical results to our disposal). 

8.2. The efFect of the IMF 

Different IMFs influence in a crucial way the injection of 
dust into the ISM. For every generation of stars, the relative 
amount of newly born massive objects able to produce dust 
via the SNa channel and the amount of low and interme- 
diate mass stars refueling the ISM with the dust produced 
through the TP-AGB phase, both depend on the IMF. This 
last one also determines the timescales of Stardust injection 
and the relative contribution by SNae and AGB stars. 

In Fig. [S] we compare the degree of dust enrichment ob- 
tained using nine different IMFs chosen among those widely 
used in literature. The reference case is always given by the 
model QVABBCBBB. All the others are obtained from this 
by varying the first parameter of the list from case A to 
case I. First of all, let us examine the effects of the IMF 
during the first evolutionary stages, when the dust enrich- 
ment is mainly due to SNae, which also supply the seeds 
and metals for the accretion process in the ISM. We begin 
with the the solar vicinity displayed in the left panel of Fig. 
[51 We note that, at varying the IMF and keeping fixed all 
the other parameters of the model, the age at which the 
dust enrichment by the accretion in the ISM becomes com- 
parable to that by SNae can vary by about ^ 0.5 Gyr. The 
time difference can be easily explained as due to the differ- 
ent percentage of massive stars exploding as SN« (and thus 
refueling the ISM) with the different IMF (see the entries 
of Table H. 

Some IMFs (like those by Kroupa, Larson SoNe and Scalo) 
predict a small number of SN« compared to others (like 
those by Larson, Kennicutt or Chabrier) that are more 
generous in the number of massive stars and hence Type 
II SNae. Therefore, in the former case a small injection of 
dust by SNa explosions and a slow accretion (fewer seeds 
to disposal) are expected. The opposite holds true with the 
latter case case favouring the formation of massive stars. 
Furthermore, with the former case the time at which the 
ISM gets dust-rich by accretion is delayed with respect to 
the other case. This is also shown in the panels of Fig. |S1 



Table 4. Mass fractions in different mass intervals predicted 
by different IMFs. All IMFs are normalized to unity over the 
mass range of validity. The following nine IMFs are considered: 
Salpeter, Larson, Kennicutt, original/old Kroupa, Chabrier, 
Arimoto, new Kroupa multi-slope power law 2002-2007, Scalo 
and, finally, Larson adapted to the Solar Neighbourhood (See 
Sect. O for more details). The upper mass limit is always 100 
Mq , whereas the lower limit is chosen according to the discussion 
made in Sect. [3] Three mass intervals are considered according 
to the different Stardust factories. All masses are in solar units. 



IMF 


M<1 ' 


1 < M < 6 ^ 


M > 6 '' 


Salpeter 


0.6075 


0.2285 


0.1640 


Larson 


0.3470 


0.3568 


0.2962 


Kennicutt 


0.4094 


0.3883 


0.2023 


Kroupa (old) 


0.5948 


0.3016 


0.1036 


Chabrier 


0.4550 


0.3517 


0.1933 


Arimoto 


0.5000 


0.1945 


0.3055 


Kroupa 2002-2007 


0.6198 


0.2830 


0.0972 


Scalo 


0.6802 


0.2339 


0.0859 


Larson SoNe 


0.5614 


0.3130 


0.1256 



^Fractional mass of stars that do not contribute to the dust 
budget of stellar origin. ^Fractional mass of stars that 
contribute to the Stardust budget via the AGB 
channel. ^Fractional mass of stars that contribute to the 
Stardust budget via the type II SNae channel. 



where the time when the ISM accretion equalizes the dust 
enrichment by SNae at decreasing the relative percentage of 
massive stars in the IMF is marked. It is also evident that 
IMFs skewed toward massive stars produce much more dust 
of stellar origin. Consequently, before accretion in the ISM 
starts driving the evolution of the dust, large differences 
brought by the IMF are possible. 

So far we have examined the solar vicinity with a rel- 
atively mild star formation efficiency. What about the in- 
nermost region of the MW characterized by a much higher 
SFR? The situation is shown in the left panel of Fig. [51 
which displays the dust enrichment due to SNae as in Fig.|5l 
Compared to the solar vicinity, we note that SNas produce 
many more seeds and metals, accretion in the ISM develops 
faster and becomes important very earlier on. However, in 
the ISM the dust production by accretion becomes more 
or less comparable to that by SNae at the same time inde- 
pendently of the IMF (see the large dots marked in the left 
panel of Fig. ^ . In the case of the solar vicinity, the cross- 
over stage hardly occurred below 1 Gyr extending up to 1.5 
Gyr, whereas now they all fall in the age range 0.6-0.7 Gyr. 
In a medium rich of seeds, dust accretion grows faster and 
the effect of the IMF somehow loses importance. 
The differences both in the amounts of dust of stellar ori- 
gin injected and the timescale of earlier enrichment in dust 
by SNae are very large. In general, the effects induced by 
variations in the IMF can be very large during the earliest 
stages of evolution. Along this line of thought, we can ex- 
pect that in high-redshift obscured galaxies with high SFR 
(easily even higher than the early SFR of the inner regions 
of the MW) some IMFs may not be able to produce the 
amounts of observed dust of stellar origin before the dust 
accretion process has become significant. The immediate 
implication of t his for p r imeva l galaxies can be easily fore- 
seen. According iDraing (|2009l ). some accretion in the ISM 
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Fig. 5. Temporal evolution of the contribution to the abundance of dust during the first 2 Gyr-3 Gyr in the SoNe at 8.5 
kpc and for different IMFs. All the contributions have been corrected for the destruction of dust. Nine IMFs have been 
considered: Salpeter (thick dotted line), Arimoto (thick dot dashed line), Chabrier (thick continuous line), Kennicutt 
(thick dashed line), Larson (thin dotted line), Kroupa-2007 (thin dot-dashed line), Kroupa original (thin continuous line), 
the Scalo IMF (thin dashed line) and, finally, the Larson IMF adapted to the SoNe (thin dashed line). See Sect. [3] for 
more details. Left panel: temporal evolution of the total contribution to the dust budget by SNae in the MW SoNe. 
The filled symbols represent the instant when the contribution by ISM accreted dust equalizes the SNae injected dust 
amount. Middle panel: time evolution of the total contribution to the dust budget by AGB stars for the MW SoNe. 
The filled symbols represent the instant when the contribution by AGB injected dust equalizes the ISM dust production 
by accretion, while the empty symbols represent the instant (if eventually it happens) when AGB dust equalizes the SNae 
injected dust. Right panel: time evolution of the total dust budget for different IMFs. 
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Fig. 6. Temporal evolution of the contribution to the abundance of dust during the first 1-2 Gyr in the inner part of the 
MW disk at 2.3 kpc from the centre and for different IMFs. All the contributions have been corrected for the destruction 
of dust. Nine IMFs have been considered as in Fig. [51 See Sect. [3] for more details. The meaning of all the symbols is the 
same as in Fig. [5] 



to explain the amount of dust in observed primeval galax- 
ies is required. The picture should be as follows: (i) before 
dust by accretion in the ISM and dust injected by stars 
become comparable, the effect of the IMF prevails and de- 
termines the amounts of dust present in the Galaxy; (ii) if 
the SFR is high, the IMF does not play an important role in 
determining the onset of the dust accretion process in the 
ISM, whereas if (iii) the SFR is low, different IMFs cause 
an important spread in the ages at which dust by accretion 
becomes important. 

We pass now to examine the role played by AGB stars. 



In the middle panels of Figs. [S] and [B] we see for the same 
regions of the left panels the contribution of the AGB stars 
to the total dust budget. The AGB stars contribute sig- 
nificantly over a longer time-scale: the filled symbols show 
the age at which dust produced by AGB stars and dust ac- 
creted in the ISM become comparable, whereas the empty 
ones show the same but for AGB stars versus SNae. In both 
cases the age is marked when and if the equality among 
the three contributions can be established. There are in- 
deed some extreme IMFs, where massive stars are favored 
with respect to the intermediate mass ones progenitors of 
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AGB stars (like the Arimoto IMF), in which during the 
first Gyrs AGB stars never reach SNse in the inner MW 
regions (Fig. [5]- middle panel). Furthermore, IMFs richer 
in intermediate mass stars produce bigger amounts of dust 
by AGB stars. This may somehow be correlated with the 
delayed appearance in a galaxy spectral energy distribution 
of the PAHs features: the delay may depend on efficiency 
of dust injection by AGB stars. 

Finally, in the right panels of Figs. ([5]) and (jH]) we show 
the evolution of the total dust content in the ISM: it is in- 
teresting to note how the differences between the various 
IMFs in the early stages are the same as those expected 
(observed) at the present age, keeping constant all other 
parameters. The IMFs producing more dust by SNas (and 
of course seeds for accretion) in the early stages are the 
same for which we get higher amounts of dust at the cur- 
rent age. The differences among the various IMFs can be 
significant. Indeed, a fast enrichment of dust during the 
early stages caused by SNae, (keeping fixed all other pa- 
rameters, the condensation coefficients in particular) goes 
together with a strong enrichment in metals. More metals 
means more atoms available for dust to grow in the ISM. 
For this reason IMFs skewed toward massive stars favor 
the accretion of dust in the ISM and lead to higher final 
contents of dust. 



8.3. The effect of the SF law 

As presented in Sect.[3l four SF laws have been considered. 
The results are shown in the four panels of Fig. [7] limited 
to the early stages of the evolution. The reference model is 
gVABBCBBB, in which the SFR is changed as indicated. 

The parameters k and v of t he SFR are chosen in agree- 
ment with the analysis made bv lPortina ri & Chiosi (199!|), 
in such a way that the sole effect of the SFR law is isolated. 
All other parameters of the model are kept fixed. At given 
k and i^, the different SF laws widely adopted for the MW 
disk produce a similar dust budget. Some differences can 
be noted (i) in the amount of dust before the onset of the 
ISM accretion, with the Schmidt and/or Wyse & Silk laws 
favouring a higher amount of dust by SNae and (ii) in the 
age at which the production of dust by accretion in the ISM 
becomes more important than that by SNs. Depending on 
the SF law, the time interval in which SNae dominate the 
total dust budget gets short or long. 

Given that the specific expression for the SF law is not 
of primary importance here (at least choosing among the 
ones we included in this study), we turn the attention to 
the efficiency of SF represented by the parameter v. In Fig. 
[S] we show the contribution to the dust budget in the usual 
three significant rings of the MW at varying i/ from i/ = 0.30 
to i^ = 0.70. As expected lower values of v imply a smaller 
rate of SNae and number of AGB stars, therefore a delay 
in the ISM dust accretion process because there are less 
seeds/metals injected by stars into the ISM from the stars, 
and finally a lower total production of dust. As expected, 
varying the SF parameters affects the system in the early 
phases of the evolution when the SF is strong. This can be 
noticed once comparing thin and thick lines in Fig. [5] 

Does the efficiency z/ affect also the dust budget at the 
present epoch? In Fig. IHlwe show the evolution of the same 
three regions of Fig. |8] up to the present age of the MW. 
As we can see, in the inner regions and the solar vicinity, 
even if there is a significant difference in the past, at the 
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Fig. 7. Temporal evolution of the dust budget in the SoNe 
during the early stages until 1.5-2.5 Gyr. Four SF laws are 
shown in the four panels: Schmidt, Talbot & Arnett, Dopita 
& Ryder and Wyse & Silk. We display: the amount of dust 
grains accreted in the ISM (continuous lines), the amount 
of dust injected by AGB stars (thin dot-dashed line), the 
total amount of SNae Stardust (dashed line) also subdivided 
into the various grain families, respectively, i.e. dotted line 
(iron-dust), silicates (continuous line), carbonaceous grains 
(dashed line), and S/Ca/N based grains (dot-dashed line). 
Left panel: Temporal evolution of the dust budget injected 
in the SoNe by the various sources and the Schmidt law. 
Central-left panel: the same as in the left panel but for 
the Talbot & Arnett law. Central- Right panel: the same 
as in the left panel but for the Dopita & Ryder SF law. 
Right panel: the same as in the left panel but for the 
Wyse & Silk SF law. 



present time the difference gets negligible, whereas in the 
outer regions the difference in the dust budgets (both total 
and partial ones) remains remarkable. This is an effect of 
the adopted SF laws that are all scaled to the current star 
formation at the solar neighbourhood and for this reason 
tend to ultimately produce the same result as the evolu- 
tion proceeds. However, we expect that if the SF law is not 
tied up to normalization or scaling factor, the adoption of 
different SF laws would have a strong impact on the whole 
evolution. The difference in the outer regions is explained 
by the long delay in onset of the ISM dust accretion pro- 
cess in the case of low efficiencies v (very low star forming 
environment): when eventually the ISM accretion process 
becomes important, there is not enough time to reach the 
dust budget produced in the high v case. This behaviour 
is also strengthened by the inward radial Sows that remove 
gas from the outermost regions. 

8.4. The effects of different models for dust accretion in the 
ISM 

In this study we consider two models of dust accretion 
in the ISM (see Se c ts. 14.11 l42l and [5] for details) : model 
A based on iDwekl ()1998[ ) where accretion is simply in- 
cluded with a general timescale depending on the destruc- 
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Fig. 8. Temporal evolution of the contribution to the dust 
budget during the first Gyrs at varying the coefficient v 
of the star formation law, namely the Dopita & Ryder SF 
law, for three regions of the MW. The thin line is i^ = 
0.30 (low efficiency), while the thick line is for v = 0.70 
(high efficiency). We show: the total amount of dust (from 
accretion in the ISM plus dust ejected by SNae and AGB 
stars - continuous line); the contribution by accretion of 
dust grain in the ISM (dotted line); the total contribution 
by SNae (dashed line); the contribution by AGB stars (dot- 
dashed line). Left panel: the results for an inner ring of 
the MW. Central panel: the same as in the left panel but 
for the SoNe. Right panel: the same as in the left panel 
but for an outer ring of the MW. 
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Fig. 9. Temporal evolution of the contribution to the dust 
budget up to the present age at varying the coefficient v 
of the star formation law from 0.30 to 0.70. The meaning 
of the symbols is the same as in Fig. [51 Left panel: the 
results for an inner ring of the MW. Central panel: the 
same as in the left panel but for the SoNe. Right panel: 
the same as in the left panel but for an outer MW ring. 
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Fig. 10. Temporal evolution of the dust production during 
the early phases of the MW evolution at varying the ac- 
cretion model used to describe dust formation in the ISM. 
Three regions have been considered as usual: an inner one 
(left panel), the SoNe (middle panel), and an outer region 
(right panel). Thick lines represent model A based upon 
iDwekl (|1998[ ) and iGalura et al (2008) , wherea s thin hues 
represent model B based upon [Zhukovska et al.l ()2008() . For 
AGB (dot-dashed line) and SNae (dashed line) we only 
show one line as the contribution is fixed. We show: the 
total amount of dust grains in the ISM (continuous lines) 
and the total amount of accreted dust in the ISM (dotted 
lines). Left panel: The results for an inner region. Central 
panel: The results for the SoNe. Right panel: The results 
for an outer region. 



tion timescale. Indeed, the accretion timescale is half of the 
destructio n ti mescale, and the same value is used for all the 
ele ments JD^ ek 1998; Ca lura et aLll2008D . Model B based 
on IZhukovska ct al., (.2008 ) in which a different timescale 
for each element is adopted and the evolution of the dust 
abundances for a number of elements, supposedly contained 
in a several species of dust grains, is followed. Furthermore, 
there is no a priori connection between the destruction and 
accretion timescales. 

In Fig.[TOl we show the evolution of the dust budget for 
the two accretion models. The amounts of dust produced 
by SNae and AGB stars and all the other parameters are 
kept fixed. In particular, the description of the dust de- 
struction process is the same in both cases: therefore we 
can examine the sole effect of accretion. In the very early 
stages there is no difference: the budget is dominated by 
SNae. The accretion process however starts to be signifi- 
cant very fast in model ,4, because it simply depends on 
the adopted timescale, while in model B it is more sensi- 
tive to the physical conditions of the environment: in gen- 
eral it tends to slow down at decreasing SFR and densities 
(i.e. passing from the innermost to the outermost regions). 
Even if there is some difference between model A and B^ 
the behaviour of the total dust budget is similar. This find- 
ing means that if we are interested in the total amount 
of dust produced, we can do it simply choosing a suitable 
timescale for dust formation. This holds everywhere but the 
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Fig. 11. Temporal evolution of the dust production in the 
MW up to the present age at varying the accretion model 
used to describe dust formation in the ISM. Left panel: 
The results for the inner ring. Central panel: The same 
but for the Solar Neighborhood. Right panel: The same 
but for the outer ring. 
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Fig. 12. Temporal evolution of the dust production in the 
MW up to the present age at varying the accretion model 
used to describe dust formation in the ISM. In this plot we 
put into elements the contribution of individual element to 
the total dust budget. Left panel: The results for an inner 
ring of the MW. Central panel: The same but for the 
SoNe. Right panel: The same but for an outer ring of the 
MW. 



outer regions where for model A we probably need a longer 
timescale of dust formation to closely agree with model B. 
This sounds reasonable because it is likely that the accre- 
tion timescale varies with the environment. In any case even 
the simple model of dust accretion with constant timescale 
is fully adequate to follow the evolution of the total dust 
budget in the early stages of evolution. 

In Fig. [TTj we extend the evolution up to the present 
time: in the inner regions and solar vicinity of the MW the 
differences between the two models are quite small, whereas 
as expected they are large in the external regions of the 
disk. In brief, in model A the accretion in the ISM starts 
very early independently from the environment and conse- 
quently it gives rise to a higher dust content compared to 
model B. 

A deeper insight of the differences brought by the mod- 
els of dust accretion is possible looking at the evolution 
of the single elements composing the dust. In Fig. [T^] we 
show the evolution of some depleted elements up to the 
present time, namely C, O, Mg, Si and Fe. As expected 
there is a strong disagreement between the two accretion 
models in all the regions. In the case of oxygen the differ- 
ence is striking: model A with a fixed timescale produces a 
lot of oxygen in dust. Since there is no description of how 
the various elements enter the different dust grains, with a 
fixed timescale the most abundant element is also the most 
abundant in dust. With a simple model we can not follow, 
for instance, the different ways in which oxygen is bound 
in silicates, or iron is bound both in iron-dust and silicates. 
Of course, only the comparison with the observational de- 
pletion factors can highlight the issue. Most likely, model B 
taking into account the physical conditions of the medium 
and following in detail the evolution of typical dust grains 
should better reproduce the observational depletion factors. 



8.5. A final note on the model parameters 

We have just discussed the main parameters entering the 
problem of the calculation of the dust enrichment of the 
ISM. However we must underline the following point: clas- 
sical chemical models are widely used in literature to follow 
the metal enrichment in galaxies of different morphology. 
They are already complicated because of the many phys- 
ical ingredients entering the problem, like the law of SF, 
the IMF, the stellar yields and the geometrical descrip- 
tion of the galaxy. When the dust is added to the problem, 
the parameter space literally blows up, because in addition 
to the classical parameters we have to consider also those 
governing the dust content. The amount of dust of stellar 
origin injected in the ISM, the way in which the dust con- 
tent increases/decreases by accretion/destruction and other 
physical processes, all these concur to extend the list of pa- 
rameters. In addition to those we have already discussed, a 
couple of them deserves some discussion. 

(1) The chemical model we are working with has a single- 
phase description of the gas content and in particular 
of the MC component in which the dust accretion takes 
place. For this reason we have included two possibili- 
ties for the MC fraction xmc- Either xmc is assumed 
to be constant and equal to the present day value in 
the solar vicinity, or it stands on the data for the MW 
and it correlates the fraction of MCs to the SFR and 
the gas density through AN Ns. In this la t ter op tion, 
that is described in detail in iPiovan et al.l (|2011c[ ). we 
are therefore assuming a variable amount of MCs. Both 
the recipes tend to the same value for the current time 
and produce similar results for the SoNe, that is for our 
target, since the differences in xmc^ most striking in 
the early stages of the evolution, are not enough to pro- 
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duce significant differences. The reason for this can be 
traced in the SFRs and densities of the SoNe that are 
never characterized by extreme values. Instead, quite 
different results between the two recipes are expected 
for the outer or inner regions of the MW, because of 
the more significant excursion of the physical variables 
(|Piovan et al.ll2011 j ). 
(2) All our models includes radial flows of matter and the 
effects of a galactic bar. All the details and typica l pa- 
rameters are taken from Portinari &: Chiosil (|l999l ) and 
iPortinari fc Chiosil ()200[ll) . Radial flows and bar have 
been taken into account to better reproduce a wide num- 
ber of properties of the MW disk, in particular to simul- 
taneously reproduce the radial gradients and the peak 
of gas observed around 4 kpc. The main motivation for 
including the radial flows and the bar is therefore the 
consistency between theory and observation as far as 
the depletion factors, radial/local abundances and gas 
masses are concerned. However, the bar influences the 
innermost properties of the MW and it is of no inter- 
est in this work. Again, the influence of the radial flows 
of m atter mainly applies to the radial properties of the 
disk (jPiovanet al.ll2011cD . If limited to the SoNe, their 
effect could be mimicked with a slightly different choice 
of the parameters. Concerning the radial flows, it must 
be flnally observed that in principle the pattern of veloc- 
ities of gas and dust could be different. For the sake of 
simplicity, we assume here that gas and dust are moving 
with the same velocities. Thanks to this gas and dust 
should have similar radial behaviour (see iPiovan et al.l 
I2011d. for aU details). 

9. Models and observations 

The final step of our study is to compare the theoretical 
results with the data in the Solar Vicinity, thus ultimately 
validating the model we have built up to describe the dust 
enrichment of the ISM. As we already discussed, the model 
has many parameters, which together with the many ob- 
servational data to match would make the search of best-flt 
solution a huge task to be accomplished. However, pursu- 
ing this strategy would not lead us to get a deeper physical 
insight of the dust formation in the MW Disk and Solar 
Vicinity in particular. Therefore, instead of looking for the 
absolute best-fit model, we are more interested in the model 
response to thoughtful choices of the parameters and in the 
comparison of model results with data on element deple- 
tion, local abundances. The following observational data 
for the MW are taken in consideration: 

(1) The depletion of the elements in the local ISM (see 
Sect. [7]) is the main check for a dust accretion model. 
This one to be physically consistent must reproduce the 
observational depletion of many elements. Obviously the 
depletion is line-of-sight dependent and we can only es- 
timate the range of plausible values (See Fig. [2]). 

(2) The local evolution of the elemental abundances in 
the Solar Vicinity. This is observationally indicated by 
many diagnostic planes such as [El/H] vs. [Fe/H] and 
[El/Fe] vs. [Fe/H] for some elements that are also in- 
volved in the dust formation process, derived from large 
samples of F and G stars. The models need to match 
those diagnostic planes. This allows us to check that not 
only the process of dust formation/injection is properly 



simulated but also that the total enrichment process, as 
observed in the different generations of stars, is realis- 
tically reproduced. 

(3) The large scale properties of the MW in the Solar 
Vicinity like the surface densities of stars and gas, rate 
of SN«, dust-to-gas ratio, and present day SFR must 
be reproduced. 

(4) The age-metallicity relation in the [Fe/H] vs. Age plane 
as observed in stars of the Solar Vicinity. 

(5) The metallicity of the Sun at the current age and of 
the proto-Sun about 4.56 Gyr ago. 

The reference model is always QVABBCBBB according 
to the parameter list contained in Table |3] and it includes 
the radial flows and the bar. The main parameters we 
are going to play with are the star formation efficiency 
V, the infall timescale r and the IMF. For this latter we 
consider those by Kroupa, Larson (however adapted to 
the Solar Vicinity in the high mass tail as discussed in 
Sect. [S]) and by Salpeter (usually taken as the reference 
case for comparison) . As already listed in Table 01 we limit 
ourselves to the IMFs that do not predict a high number of 
SNa8. IMFs of this type are perhaps more suited to other 
dust-rich environment s , like starbu r st gala xies or ellipticals 
(IValiante et al.l 120091: iGaU et al.l l2011allbl : iPioino et all 
1201 ll ). The MW disk and the Solar Vicinity in particular 
seem to req u ire IMFs somewha t poor in massive star s 
(iDwekl 119981 : ICahira et al.l 120081: IZhukovska et"!!! l2008l^ . 
Four infall timescales are considered, that is r = 1, r = 3, 
r = 6 and t = 9, from the shortest to the longest, and four 
values of the star formation efficiency, v = 0.3, v = 0.7, 
V = 1.1 and ly = 1.5, from the lowest to the highest 
efficiencies. 

The results of these models are presented in the vari- 
ous panels of Fig. 1131 In general, moving from the upper 
right to the bottom left corner of the figure, the depletion 
of the elements gets easier to obtain. We find that: (i) in 
the case with the Kroupa IMF the observational range of 
abundances is easily reproduced over ample ranges for the 
remaining two parameters (v and r); (ii) models with the 
other two IMFs lead good results only if the infall timescale 
T and the efficiency v are properly chosen. About the de- 
pletion of the individual elements we note what follows: 

Carbon: carbon depletion is simulated assuming that 
the fraction of carbon hidden in the CO molecules amounts 
to 0.3. This parameter plays a key role as shown in Fig. [31 
At varying ^co it is possible to allow for less or more dust 
embedded in C grains, whose accretion depends on the C 
atoms free to accrete. Carbon depletion is, in any case, well 
reproduced for the most common sets of the parameters. 

Nitrogen: to reproduce the small depletion of nitro- 
gen we find that the simple choice of the longest timescale 
between the oxygen and carbon, and the use of a mean ni- 
trogen dust grain were not enough to keep low the amount 
of nitrogen condensed into dust. To this aim, we introduce 
a multiplicative factor Nx for the accretion time scale t^ j^ 
that enters Eqn. (|7T|) and tune it so that the low observa- 
tional depletion is reproduced. The results shown in Figs. 1131 
and [T4l include the scaling factor Nx- We get that Nx > 10 
is required to fetch the small accretion of nitrogen. 

Magnesium and Silicon: magnesium is embedded 
into dust thank to the presence of olivines/pyroxenes and in 
our best simulations it is found to be easily very depleted. 
The preference goes toward the maximum allowed depletion 
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Fig. 13. Depletion of C, N, O, Mg, Si, S, Ca and Fe in the ISM as observed in the SoNe. The observations are compared 
with the models at varying three important parameters, namely the IMF, the SF efficiency v^ and the mass accretion time 
scale T. For the IMF we consider : the recent Kroupa IMF (left panels), the classical Salpeter one (central panels) and 
the Larson IMF adapted to the SoNe (right panels). These IMFs are all described in Sect. [21) • Four cases are considered 
for I/, that is ^ = 0.30 (empty circles), y = 0.70 (filled circles), v — 1.10 (empty squares) and v = 1.50 (filled squares). 
Finally, four values are used for the accretion time scale, namely r = l,r = 3,r = 6 and r = 9 Gyr , from the top line 
to the bottom line of each panel. 



rather than the minimum one. The opposite happens for sil- 
icon which is more often found close to the minimum deple- 
tion limit (see Fig.lT^. This common behavior for Mg and 
Si can be explained with the u nder-production of Mg char- 



acterizing the SNae yie lds by Portinari et al. (I1998D that 



are ultimately based on Wooslev fc Weaved (I1995D: This is 
a lon g known problem ()Timmes et al.|[l995t [Francois et al.l 



l2004 l that becomes crucial in our case, because Mg is one 
of the possible key elements driving the accretion process 
of the silicates. If the amount of Mg atoms available to 
form dust is under- abundant, this will set up an upper limit 
to the amount of olivines/pyroxenes that can be formed. 
In this study, first we calculate models with the original 



yields by iPortinari et al.l ()1998[ n (filled and empty circles 
and squares in Fig. [T3l) . This choice leads to a dust mix- 
ture where Mg (already under-abundant) is usually more 
depleted than Si. The case is w -ell ill ustrated in Fig. [T^ 
Accordin g to IZhukovska et a .1 (I2008D . the yields of Mg 
based on lWooslev fc Weaver! ()1995[ ) could be suited to re- 
produce the Milky Way once scaled by a factor of 2-3 in 
order to obtain Mg/Si/Fe ratios in agreement with the ob- 
servational ones for the MW, thus leading to a dust mixture 
closer to reality. Keeping in mind this suggestion, we cal- 
culate the same models of Fig. [13] slightly modifying the 
yields of Mg for Type II SN g;, in practice changing the 
tabulations of IPortinari et all (1998) according to the sug- 



^ The latest version of the yields is used in which the under- 
abundance of Mg is partially corrected (Portinari, 2011 private 
communication) . 
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gestion bv lZhukovska et aTJ (|2008[ ). These new calculations 
allow us to test the model response to variations of the Mg 
yields. The results are shown in Fig. [T3| where we note that 
the agreement with the observational data is better for all 
the IMFs, in particular with the Kroupa IMF the results are 
very good. Furthermore, the whole range of observational 
depletions is now better covered by varying the star forma- 
tion efficiency than with the original yields of Mg: changing 
V from 0.3 to 0.7 all the range of observational values is ob- 
tained. However, as in the previous case of Fig. |T21 the 
highest star formation rates enrich too much the ISM in 
metals and the process of dust formation is hardly able to 
deplete the ISM. The models clarify that the process of dust 
formation crucially depends on the the amounts of metals 
injected by the stars via mass loss and/or SNae. They also 
make clear that the set of abundances at the base of the 



theoretical models for metal enrichment and dust produc- 
tion of the ISM must be strictly identical, otherwise there 
would be no consistency between the two descriptions (con- 
sistency is of course always secured in reality) . 

Iron: in normal circumstances iron is highly depleted. 
In most of the models in Fig. |T3| the best we can obtain at 
varying the parameters is to reach the upper limit of the 
observational depletion range, i.e. —2 < [Fegas/-ff] < — 1- 
Even if two possible processes for the formation of dust 
(iron dust and silicates) are included, the mechanism of 
iron formation in cold regions of our models is not able to 
reach such lower values. Probably a more complex model, 
for instance with a spectrum of MCs with different life- 
times and/or some processes of accretion in other parts of 
the ISM different from MCs, is needed to reach the severe 
—2 depletion limit. With the modified Mg abundances, the 
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Table 5. Comparison of the properties of the MW in the SoNe with model results^ and with the IZhukovska et alJ (|2008l ') model 
from which the observational data is taken: 



Observable 



Observed ZGTOS^ 



This work 



Reference 



iHolmb erg fc Flvnnl (|20"oi r 
Dickcv (1993) 
IDam c (1993) 

lOUing fc Merrifieldi (pOOH ) 

iGilmore et alJ fioM) 

iMera ct al. (1991) 

Rana (1991) 

Ir^mmann et al. (1994) 

iTammann et al.l (|1994l ) 

Braun fc Thilker (2004) 



Total surface density a {j-Qytc) [Mqpc 
ISM surface density ctai (^©,^0) [Mqpc" 



Gas fraction (Tai (r©, ic) /o" (r0,fG) 
Surface density of visible stars cr* {rQ,tG) [Mqpc" 
Surface density of stellar remnants [Mqpc"^] 
Star Formation Rate [M0pc~^Gyr~^] 
SNffi II rate [pc-^Gyr" 
SNffi la rate [pc~^Gyr~ 
Infall Rate fMepc-^Gyr^^^l 



50- 62 



7 



13 

8 

- 14 

- 0.2 
40 



13- 
0.05 
30- 
2- 4 
3.5 - 5 
0.009-0.0326 
0.0015-0.0109 
0.5 - 5 



56 

9.7 

9.7 

9.7 

0.17 

38.6 

7.7 

3.1 
0.016 
0.0024 

1.45 



52 



10- 

10- 

10- 

0.2- 

23 

6- 

1.6 

0.01 

0.0035 ■ 

0.25- 



19.5 
19.5 
19.5 
0.39 

- 35 
10 
-4 

- 0.04 

- 0.0078 

- 1.80 



^The range of values for the entries of column (4) refers to the models discussed in this section, for the case with no correctio n to 
Mg ab undance. Of all the models shown in Fig. 1131 we select only the twelve models presented in Fig. 1151 IZhukovska et ahl 
l(200l . 



results improve a bit due to the higher amount of silicates 
and hidden iron, but still we are not able to fully span the 
interval from —1 to —2. 

Calcium: some of the considerations we made for the 
iron apply to calcium, for which the observed depletion is 
even more severe than that of iron. For calcium, there is no 
overlap at all with the observational interval. As described 
in Sect. 15.61 we introduce a multiplicative factor Cax in the 
accretion time scale t^^ ^^ and try to calibrate it against 
the observational values. Even using Cax ^ 0.05 it is not 
possible to describe the strong depletion of Ca. The results 
presented in Figs. [T3] and [U] are for a small scaling factor 
~0.1. 

Sulfur: sulfur can be very depleted along some lines of 
sight, much less along other. Since a typical sulfur dust 
grain is not available to simulate the accretion process, 
we used a simple prescription for the average accretion 
with sulfur atoms accreting on themselves (see Sect. 15.71 
for more details). A multiplicative factor Sx has been in- 
troduced to calibrate the timescale Tg" , the longest between 
the timescales of the refractory elements and the accretion 
of sulfur atoms on themselves. The results displayed in Figs. 
[T51 and [m show the sulfur depletion with (filled and empty 
hexagons) and without (filled and empty circles) the ef- 
fect of Sx on the accretion time scale Tg" . With Sx going 
from 1 (no correction - circles) to 0.6-0.4 (small correction - 
hexagons) the theory fairly agrees with observational data. 

Oxygen: this element does not participate as a key- 
element in any accretion process of dust formation. However 
it takes part to the formation of the silicates and it is an 
ingredient of the dust yields from stars of different masses. 
Because of its high abundance, the depletion is small, even 
if lots of oxygen atoms are contained in olivines/pyroxenes. 
In Figs. [13] and [T3] we can see that the theoretical predic- 
tions well reproduce the observational data, in particular 
when Mg is corrected for under- abundance in the yields. 
The agreement simply follows from using the correct Mg/Si 
ratio. 

Abundances and Depletions in the Solar Vicinity. To check 
the internal consistency of the depletion models, we must 
secure that the theoretical evolution of the abundances 



of the various elements in the Solar Neighborhood is able 
to reproduce the observational data. It would be a point 
of strong contradiction if the models can reproduce the 
elemental abundances in the dust, but fail to reproduce 
the pattern of abundances in the gas and stars. In Fig. 
115! we show the evolution of the metallicity Z and iron 
abundance [Fe/H] for a wide selection of models, taking 
into account more or less the combinations of IMF, v and 
T that best reproduce the depletion measured in the SoNe. 
In such a case, one should discard the models with r = 1 
that hardly fit the data for depletion in the SoNe. Most of 
the remaining models fairly agree with the observational 
data with the exception of the case with the Kroupa IMF, 
slow infall (r = 9), and low SFR [v = 0.3) that is not able 
to grow fast enough in metals. In Table [5] we present the 
check of the minimum/maximum values obtained in our 
models vs. the observations in the SoNe of various physical 
quantities of interest. The same simulations presented in 
Fig. [15] are compared with both the results obtained by 
iZhukovska et al. (2008 ) and the observational data. The 
general agreement between our results and the observations 
is good and it allows to conclude that we are employing a 
correct modelling of the SoNe. Once that the more general 
quantities, like the mass of stars, gas, the SNa rates and 
the total surface mass are satisfactorily reproduced we can 
proceed to examine the evolution of the abundances of the 
single elements in the SoNe and how they match up with 
the local data. 



In the panels of Fig. [16] we show the time evolution in 
the Solar Neighborhood of the abundances of six elements 
heavily involved in dust formation, namely C, O, Mg, Si, S 
and Ca, in the diagnostic planes [A/Fe] vs. [Fe/H]. The ob- 
servational data refers to F and G stars (see Fig.[TB]for more 
details on the data and associated legend). Twelve models 
are displayed at varying the IMF (Larson SoNe, Salpeter 
and Kroupa), efficiency of star formation {v = 0.3 and 0.7), 
and the infall timescale r (3, 6 and 9 Gyr). Not all the cases 
are shown for t — Z and r = 6 but only the most interesting 
ones as far as the degree of depletion is concerned (see Figs. 
[T5I and [n)l . In general, the agreement is good for t = 3 and 
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oretical abundances fairly agree and consequently dust for- 
mation stands on a realistic description of the evolution of 
the elemental abundances in the Solar Neighborhood. 

10. Discussion and conclusions 



[Fe/H] 

Fig. 17. Evolution of the Mg abundance in the Solar 
Neighborhood as measured in a sample of F and G stars. 
The yields have been slightly corrected for the Mg under- 
abundance. This corrected abundance is the same as in Fig. 



M 



T = 6, while for r = 9, during the first 2 Gyr (in practice 
until [Fe/H] ~ —1.5) the enrichment in metals is too slow 
with respect to the observations. The only exception is car- 
bon which always agrees with the observations. It must be 
pointed out that we are using infall models with a single 
timescale. This simplified picture leads to models that in 
the [A/Fe] vs. [Fe/H] planes cannot reach regions of very 
low or very high [Fe/H]. In reality, during the lifetime of 
the MW different time scales could be involved in the evo- 
lution of the describe the Solar Neighborhood, e.g. a fast 
early enrichm ent followed by a much slower one . Double 
infall models (|Chiosil[l980t IChiappini et al.l[l997f) or more 
complicated scenarios are however beyond the aims of this 
study. Finally, it is worth noticing; (i) the under-abundance 
of Mg in the top right panel; (ii) the marginal agreement 
of the models with sulfur evolution. 

In Fig. [T7]we show the evolution of the Mg abundance in 
the SoNe once corrected for the under-abundance problem. 
The same models already presented in Fig. [TC] are shown, 
but for the correction of the Magnesium abundance. With 
this correction not only we have a better pattern of deple- 
tions (See Fig. [H]), but also a better reproduction of the 
local data. 

To enforce the argument we examine the evolution of 
the most important elements involved in the dust formation 
process looking at the [A/Fe] vs. [Fe/H] diagrams relative 
to the SoNe. The most striking points of disagreement with 
the observational data are the under-abundance of Mg and 
the evolution of Sulph ur. For the sake of com parison, we 
also show the models bv lZhtikovska et aL (|2008|) , calcu lated 
after revising t he yie lds by IWooslev fc Weaven (|1995[) and 
iNomoto et al.l (|2006f ). In both cases the yields of Sulphur 
do not lead to a good fit of the observations. The yields 
bv lNomoto et al.l (120061 ) give indeed the worst model. Our 
yields produce similar results in the sense that they tend to 
keep the abundance of Sulphur too high. Fortunately, this 
has no effect on the overall production of dust as Sulphur 
drives its own accretion efficiency (See Sect. 15.71) and does 
not affect other channels of dust production. Apart from 
these minor uncertainties, the observational data and the- 



In the first paper of this series of three (jPiovan et al.l 
l2011a[ ). we presented and discussed the prescriptions cur- 
rently in use to describe the type and amounts of dust in- 
jected into the ISM by AGB stars and SNa explosions. The 
condensation coefficients of the dust have tested and ana- 
lyzed referring to a suitable chemical model for the MW 
Disk and the Solar Neighbourhood in particular. This ref- 
erence model and its physical ingredients are described in 
great detail in this paper. In particular, we focused on the 
mechanism of dust accretion in the ISM, dust destruction 
by various processes, the main parameters of the chemi- 
cal model, namely IMF, star formation efficiency and infall 
timescale. The main conclusions can be summarized as fol- 
lows: 

— The CO molecules influence the formation of C-based 
dust, thus introducing into the estimates an unavoidable 
uncertainty. The higher the amount of carbon embed- 
ded into CO, the slower is the accretion process and 
the smaller the amount of carbonaceous grains that are 
formed. 

— During most of the MW evolutionary history, the main 
process enriching the ISM in dust is the accretion in the 
cold regions. Only in the very early stages, SNa3 dom- 
inate, and the duration of this phase tends to shorten 
in regions of high star formation and fast enrichment in 
metals (innermost regions of the MW). The opposite for 
the regions of low star formation (the outskirts of the 
MW), where the accretion in the ISM becomes signifi- 
cant much later so that for many Gyr SNae govern the 
total dust budget. In this case, AGB stars play an im- 
portant role, because there is enough time for them to 
significantly contribute to the dust budget, without be- 
ing overwhelmed by the SN« (earlier phases) or the ISM 
accretion (later phases) . The time interval during which 
SNae are the main dust producer can slightly change at 
varying the upper mass limit of stars undergoing the 
AGB phase. 

Our conclusions for the MW Disk could be extended to 
galaxies characterized by continuous star formation on 
the notion that the outer regions of the MW Disk might 
correspond to low star forming galaxies and the inner 
ones to high star form hig objects. Very high values of 
SFR (See for example IcSl et al.ll2011al lhl. for detailed 
simulations of starburst galaxies and QSOs), which are 
not reached in the MW, not even in the early phases 
of the evolution in the inner regions, somehow elude 
this simple scheme. In any case, star-dust dominates 
the mild SF environments for a long period of time. 
This time scale tends to decrease in environments with 
IMFs skewed toward the mass interval in which the star- 
dust/metals factories (AGB stars and SNs) are impor- 
tant, thus helping the onset of the accretion phase in 
the ISM. 

— In the high SFR/high metallicity regions AGB stars 
mainly produce silicates, whereas in the low SFR/low 
metallicity ones, carbon stars can contribute signifi- 
cantly to the C-based dust. 
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Fig. 18. Evolution of the elemental abundances in the Solar Neighborhood as measured in a sample of F and G stars. 
The evolution of [A/Fe] vs. [Fe/H] i s shown for six e l ement s of interest, that is C, O, Mg, Si, S and Ca. F or the sake 
of comparison we plot the results by IZhukovska et al.l (|2008[ ) for the SNa yields by IWooslev fc Weaven (|1995[ ) where Mg 
and Fe have been co rrected to get a better agreement with the observations (the thick dotted line), and the yields by 
iNomoto et al.l (J2006I ). the thin dotted line. The meaning of the symbols for the observational data is the same as in Fig. 
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The IMF plays a fundamental role because it controls 
the relative amounts of low, intermediate and massive 
stars (see the entries of Tabic 2]) . We found that IMFs 
skewed toward massive stars are not suitable to repro- 
duce the properties of dust in the MW, the SoNe in 
particular, (see als o lDweklll998l: IZhukovska et al.ll2008t 
ICalura et al.ll2008^. whereas it seems that for star-burst 



galax ies (iGah et ah 2011a), ellipticals (JPipino et all 



galax ies iHjaii et ai.i i/uiia ), empticais iinpmo et ai.i 
[Ipn') andQSOs (iGall et al.ll20TTb[|Valiante et al.ll201lD 



IMFs more biased toward higher mass stars are re- 
quired, because otherwise it would be more difhcult to 
reproduce the amount of dust observed already on site 
in very high-z objects. The alternative is to introduce 
very high condensation efficiencies in SNae. This pos- 
sibility is somehow supported by the recent FIR/sub- 
mm o bservations on the SN 1987A bv iMatsuura et aP 
(|2011[ ). In such a case normal IMFs in high-z galaxies 
and QSOs cannot be excluded. The use of different IMFs 
for different galaxies in order to reproduce dust proper- 
ties is a subject of vivid debate, and tightly related to 
the wider question about the universality of the IMF. 
The debate seems to favour an IMF sensitive to the 
initial conditions of star formation (see for instance the 
recent lElmegreenll2009t iBastian et ani2010t iMvers et all 



l201ll:lGunawardhanaet"aLll201ll:lK"roupall201lL to men- 
tio n a few). Th is possibility was suggested long ago 
by IChiosi et al.l ([1998) to solve the apparent contra- 
diction between the spectral and chemical properties 
of early type galaxies. S imilar conclusions are reached 
bv lValiante et al.1 ()2011[ ). trying to match the observed 
properties of the QSOs SDSS J1148-h5251, where a top- 
heavy IMF allow a more coherent match between differ- 
ent observations. Finally, when the SFR is high (inner 
regions), the effect of the IMF on the dust budget tends 
to disappear at the current time, because the accretion 
in the ISM becomes dominant early on in the evolution, 
while if the SFR is low we can see the effect of different 
IMFs spreading until the current time. 
Lower efficiencies v of the SF, in the range adopted for 
the SoNe, correspond to a slower onset of the accretion 
in the ISM and a lower final dust budget. 
We tested different descriptions of the process of dust 
accretion in the ISM. A simple approach adopts an av- 
erage time-scale of accretion to estimate the total dust 
budget at least in the early stages and for normal star 
forming environments. More complicated descriptions 
are require to follow the evolution of the abundance of 
single elements in dust or low star forming regions. 
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— The range of depletions observed for most of the ele- 
ments, that is C, N, O, Mg, Si and S, is nicely repro- 
duced by reasonable combinations of the parameters. 
IMFs with the classical slope in the high mass range, 
more easily fit the data for depletion. The only excep- 
tions are: (i) iron, for which probably a more complex 
mechanism of accretion in t he ISM than the simpl e one 
on cold regions is required ()Zhukovska et al.ll2008D : (ii) 
calcium whose extreme depletion can not be reproduced 
by our model. 

— The ratio between the abundances of Mg and Si is cru- 
cial in the formation of silicates and it is important to 
tune it according to the ratio expected from models. 
We had to correct the slightly under-abundance of the 
Mg in our gaseous yields, to obtain a better agreement 
with the observed depletion for the refractory elements 
involved into the silicates formation. 

To conclude, the classical chemical models nicely 
reproduce the observed depletion and properties of the 
SoNe. The theoretical ingredients, like condensation 
coefficients and accretion models, behave in a satisfactory 
way. What can be improved? Clearly the major weak 
point is the one-phase description of the ISM. At least a 
two-phases ISM is required. Finally, work is in progress to 
introduce a multi-phase description of the ISM with dust 
in N-Body-TSPH simulations of galaxies. 
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